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Abstract 

We have calculated the leading-twist next-to-leading order (NLO), i.e., correction 
to the light-cone sum rules prediction for the electromagnetic form factors of the nucleon. 
We have used the loffe nucleon interpolation current and worked in Mtv = approxi- 
mation, with Mtv being the mass of the nucleon. In this approximation, only the Pauli 
form factor F2 receives a correction and the calculated correction is quite sizable (cca 
60%). The numerical results for the proton form factors show the improved agreement 
with the experimental data. We also discuss the problems encountered when going away 
from Mtv = approximation at NLO, as well as, gauge invariance of the perturbative 
results. This work presents the first step towards the NLO accuracy in the light-cone 
sum rules for baryon form factors. 
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1 Introduction 

Exclusive processes offer challenging ground for quantum chromodynamics (QCD) and, 
especially the ones involving hadron form factors, provide us with valuable insight in 
the internal structure of composite particles. The simplest probe is the photon and thus 
obtained electromagnetic form factors characterise hadron's spatial charge and current 
distributions. 

The framework for analyzing exclusive processes at large-momentum transfer within 
the context of perturbative QCD (pQCD) has been developed in the late seventies 
[I1I21E1I11E1EIEIE1E]- it was demonstrated, to all orders in perturbation theory, that 
exclusive amplitudes involving large-momentum transfer, i.e., so-called, hard-scattering 
amplitudes, factorize into a convolution of a process-independent and perturbatively in- 
calculable soft part, i.e., distribution amplitude (one for each hadron involved in the 
amplitude), with a pro cess- dependent and perturbatively calculable elementary hard- 
scattering amplitude. In the leading-twist approximation of the standard hard-scattering 
approach, hadron is regarded as consisting only of valence Fock states, and transverse 
quark momenta are neglected (collinear approximation) as well as quark masses. In 
this picture each hard gluon exchange brings factor a^/vr, while higher-twist effects are 
suppressed by l/Q"^, with Q"^ being the characteristic large scale of the process (i.e., 
in the case of electromagnetic form factors, that is the virtuality of the photon). Al- 
though this pQCD approach undoubtedly represents an adequate and efficient tool for 
analyzing exclusive processes at very large momentum transfer, its applicability at ex- 
perimentally accessible momentum transfers has been long debated and attracted much 
attention. Even in a moderate energy region (a few GeV) soft contributions (resulting 
from the competing, so-called, Feynman mechanism) can still be substantial, although 
the estimation of their size is model dependent. Recently, the concept of generalized 
parton distributions (GPDs) (TUl HU [T2] has been introduced to describe the soft part in 
various exclusive processes (hke deeply-virtual Compton scattering, deeply-virtual elec- 
troproduction of mesons . . . ) and make a connection between inclusive and exclusive 
processes and corresponding characteristic quantities like parton distribution functions 
(PDFs) and form factors. Although more general, that approach (for details, see reviews 
[131 [E]) is basically similar to previously described pure pQCD approach, only at, for 
example, leading-twist the hard-scattering part does not involve all Fock state partons 
but instead one uses the so-called "hand-bag" picture. 
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The QCD sum rule approach [TSl [IS] apphed to the pion form factor supports the 
conclusion that the soft contributions are dominant at moderate momentum transfers up 
to ~ 2— 3 GeV^ [T71|TH]. The application of the method at higher faces the problem 
of an ill-behaved series in Q^/M|, where is Borel parameter. Moreover, for nucleon 
form factors the QCD sum rule approach only works in the region of small momentum 
transfers < 1 GeV^ |T9| [20] . One can find in the literature many approaches and 
attempts to circumvent these problems. 

The light-cone sum rule (LCSR) approach [211 [22], adopted also in this work, can 
be regarded as a successful technique which combines sum rule principles with pure 
perturbative QCD approach advocated earlier. The domain of validity extends above few 
GeV^. In LCSR approach the "soft" contributions to the form factors are calculated in 
terms of the same DAs that enter the pQCD calculation and there is no double counting. 
Hence, although LCSRs do involve a certain model dependence, the important advantage 
of this approach lies in the fact that it is fully consistent with pQCD. In the last years 
it has been widely applied to mesons, see [23], [21] for reviews. Moreover in Refs. [25l 126] 
LCSR approach was introduced for the description of nucleon DAs and nucleon form 
factors, and further analysis followed in Refs. (23, [2H] and [221 [30]- The weak decay 
Afe — »• pli^i was considered in [21] and the A^7A transition form factor was worked out in 
[S]. 

As explained in detail in, for example, Ref. [2H], the basic object of LCSR approach 
to, say, nucleon form factors, is a correlation function expressed in terms of the matrix 
element of the time ordered product of the current of interest (in our case, electromagnetic 
current) and a suitable nucleon interpolation current. The matrix element is taken 
not between the vacuum states but between vacuum and nucleon state \N{P)), which 
represents the second nucleon in the process. When both the virtuality of the photon = 
—Q^ and the momentum flowing through the nucleon interpolation current vertex P'^ = 
{P — qY are large and negative, one can use the operator product expansion (OPE) on the 
light cone, i.e., one employs pQCD to evaluate the Wilson coefficients, while the matrix 
elements of the relevant composite operators correspond to the appropriate moments of 
the nucleon DAs. This procedure is quite analogous to the determination of the hard- 
scattering amplitudes in the pure pQCD approach discussed above. Furthermore, in 
order to access the nucleon form factors, one then makes use of the dispersion relation 
in P'^ and define the nucleon form factor contribution through the, so-called, interval 
of duality or continuum threshold Sq. The usual application of Borel transformation 
facilitates further the calculation. 

There are some important features of this approach to be stressed. The leading-order 
(LO) contribution to the form factor is a purely soft contribution which is represented as a 
sum of terms ordered by twist0of the operators, i.e., nucleon DAs. Moreover, in contrast 
to the pQCD hard-scattering approach, the contributions of higher-twist DAs are not 
suppressed by but by powers of P'^, i.e., by powers of Borel parameter ^1 — 2 
GeV^. Thus their role is more pronounced. Furthermore, the LCSR expansion contains 
terms generating also the asymptotic pQCD contributions. For the pion form factor the 
hard-scattering contributions appear at order a^, and in Refs. [331131] it was explicitly 
demonstrated that they are correctly reproduced. For the nucleon form factors they 
appear at order aj.. 



^In this work under twist we assume light-cone and not geometric twist. 
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Let us now turn to the next-to-leading (NLO) contributions. It is well known that, un- 
like in QED, the leading-order predictions in pQCD do not have such predictive power, 
and that higher-order corrections are important. Still, although the LO predictions 
within the hard-scattering approach (as well as, GPD based approach) have been ob- 
tained for many exclusive processes, only a few processes have been analyzed at the NLO 
- see the detailed account in, for example, Ref. [35], and additionally Refs. [36l |371 [38] . 
Similarly, as it was stressed in, for example, [2S] the LO LCSRs may not be sufficiently 
accurate. The radiative gluon corrections to LCSRs were calculated for number of pro- 
cesses involving mesons, i.e., pion form factor [33], pion transition form factor [39j, the 
decay 5 ^ vrez/ [iQl HI], 5 to vr (iT, r]) form factor [121 Hal El US], and 5 to light- 
vector meson (p, u, K*, 0) form factors [IHIHT]. The radiative corrections to nucleon 
form factors have not been evaluated either in hard-scattering picture nor in the LCSR 
approach. 

In this work we took a task of calculating the NLO corrections to LCSRs for nucleon 
form factors. We follow closely Ref. l28j and extend the formalism to NLO calculation. 
Even at LO the LCSR formalism for baryons is considerably more cumbersome than for 
mesons. As we shall show, at NLO this is even more pronounced. For example, while 
for mesons at next-to-leading twist, i.e., twist-3, the use of asymptotic DAs along with 
Windzura-Wilczek approximation ensured the cancellation of the coUinear singularities 
without explicitly knowing evolution kernels (see, for example, jlS]), the nature of nu- 
cleon form factor calculation, nucleon DAs and corresponding asymptotic forms is quite 
different and does not enable similar simple cancellations. Actually, as we shall show, at 
the moment, only the approximation M^r = 0, with M^v being the nucleon mass, allows 
fully consistent NLO calculation. Hence, in this work we calculate NLO corrections to 
leading twist, that is twist-3 in nucleon case, in Mn = approximation, and outline the 
problems encountered when going away from this approximation. We furthermore ana- 
lyze the numerical importance of such NLO corrections applying them to the assessment 
of proton form factors. 

The paper is organized as follows. In Sec. [2] we give necessary definitions, introduce 
LCSR formalism and explain the preliminaries. The general LO formulas and detailed 
calculation of twist-3 and twist-4 contribution to the correlation function are presented in 
Sec. 131 The gauge invariance of these LO results is also discussed. The calculation of the 
NLO corrections to the correlation function is explained in detail in Sec. HI Section is 
devoted to developing needed LCSRs, and numerical results are presented and discussed 
in Sec. El We summarize and conclude in Sec. [71 There are five appendices devoted to 
some more technical details or summary of analytical results. In App. |A] the Feynman 
rules derived for the leading-twist calculation and employed in LO and NLO calculation 
are listed. The discussion of 75 ambiguity relevant for our NLO calculation is given in 
App. [HI Appendix [O is devoted to the summary of the analytical NLO results used in 
numerical calculations. The imaginary parts of selected functions needed in evaluating 
LCSRs are derived in App. [Dl while the selected LO twist-3 and twist-4 contributions 
and corresponding LCSRs are reanalyzed and listed in App. [El 
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2 Definitions and preliminaries 



2.1 Nucleon electromagnetic form factors 

The nucleon electromagnetic form factors are defined through the matrix element of the 
electromagnetic current by 



[N{P')\f-mN{P)) = N{P') 



N{P), (2.1) 



where Fi and F2 are Dirac and Pauli electromagnetic form factors, respectively. Here, 
P' = P — q is the outgoing nucleon momentum, while P and q are incoming nucleon 
and photon momenta, respectively. Furthermore, = —Q^ with > (for spacelike 
regime we are interested in, while the sign changes in the timelike region), and for on-shell 
nucleons = P'^ = M%. 

The Sachs form factors, i.e., electric and magnetic form factors Ge and Gm-, are 
related to Fi and F2 by 



Gm(Q') = Fi(Q') + F2(Q'). (2.2) 



At = they are normalized to proton and neutron electric charges and anomalous 
magnetic moments: 

G|(0)=Ff(0) = l, 

Gg(0)=K(0) = 0, (2.3) 



and 



respectively. 



GUO)=fip = 2.79, Ff(0) = /tp=1.79, 

GUO) = /in = -1.91 , F^{0) = K^ = -1.91 , ^ ■ ' 



2.2 Correlation function 

Correlation function, the basic object used in LCSR approach (see Fig. [T]for a schematic 
representation), is defined by 

T,{P,q) = zld'xe'^^ {0\T[r^{0)j,{x)]\N{P)) , (2.5) 

with P and q being the nucleon and photon momentum, respectively. Here is the 
electromagnetic current 

j^{x) = euu{x)jf,u{x) + edd{x)-ff,d{x) , (2.6) 
and in this work we choose the interpolating nucleon current of the form 

n{x) = e"^' [u\x) Tiu\x)\ T2 d\x) , (2.7) 
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Figure 1: Schematic structure of the hght-cone sum rule for nucleon form factors (Ref. 
[28]). The corresponding correlator function is given in f l2.5p . and P' = P — q. 

where a, 6, c are colour indices. Here the generic expression (12. 7p actually corresponds to 
the proton current, while the neutron current is obtained by replacing u ^ —d. Specially, 
for loffe current |1H] 

ri = C7A, r2 = 757\ (2.8) 

with C being the charge conjugation matrix. Other choices for ri{x) were discussed in, 
for example, Ref. [28], and loffe current was singled out as a most promising candidate 
for reliable sum rules. 

2.3 Nucleon matrix element of the three-quark operator 

Nucleon distribution amplitudes refer to nucleon-to- vacuum matrix elements of nonlocal 
operators built of quark and gluon fields at light-like separations, = 0. 
We are interested in the three-quark matrix element 

4 (0| e'^'"'u'^{aix)[aix, aox]a'au''/i{a2x)[a2X, aox]b'bd''^{a3x)[a3X, aox]c'c \N{P, A)) (2.9) 

where ^ 

[x, y] = Vexp[ig I dt{x - y)^A^{tx - (1 - t)y)] (2.10) 
Jo 

is a gauge factor which we will suppress in the following, while a,/3,7 are Dirac indices 
and P^ = Mff with M^y being the mass of the nucleon, i.e., proton mass for the three- 
quark matrix element considered above. Neutron case proceeds equivalently {u ^ d). 
Actually, in this work all results are expressed in terms of proton quantities but the 
neutron case is easy to derive from these (in the end results, it comes basically to the 
replacement e^). 

For the general Lorentz decomposition we introduce a convenient shorthand notation 

m 

4 (0| e''''K{a,x)ul{a2x)d''^{asx) |iV(P, A)) = ^^"^(M, P ■ x) , (2.11) 

i 

where JF^*) g {81,82, "Pi, "^2, Vi, . . . , Ve, .4,1, . . . , Aq, 7i, . . . , Tg} are invariant functions of 
P ■ X, while X^l and Ky^*'' are Dirac structures which can be read from Ref. [28] . 
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The structures Y^^^ contain nucleon spinor A^^ and we note that 

^ -\ -X«for^«G{5„P„A.} • ^2-^2^ 

The functions J-'^^^ are not of a definite twist and they satisfy 

u)f TD TD TD \ _ f ^^\ci2P ■ X, aiP ■ X, a^iP ■ x) for J'^*) G {Vj, Tj} 
J- [^a^P-x^a^P-x^a^P-x]-^ -T^^ia^P ■ x,a,P ■ x,a,P ■ x) ioi :F^^ e {S,,V,,A,} ' 

(2.13) 

Additional, 0{x'^) terms can be added to (12. lip but will not be explicitly considered here 
(see Ref. [28] instead). 

2.3.1 Light-cone kinematics 

For the twist classification it is convenient to go to infinite momentum frame and intro- 
duce a light-like vector by the condition 

g-2 = 0, ^2 = (2.14) 

as well as the second light-like vector 

I Ml ^ , , 

P, = P,-^z,j^^, p' = 0, (2.15) 

so that P — > p if the nucleon mass can be neglected, — > 0. The projector onto the 
directions orthogonal to p and z, is given by 

= 9f^u - —{p^^Zu + PuZ^). (2.16) 
pz 

In turn, denotes the generic component of orthogonal to z and p, and thus the 
photon momentum can be written as 

P ■ q 

qi, = q±i, + z^- — , (2.17) 

P ■ z 

where the use has been made oi p ■ q = P ■ q and p ■ z = P ■ z. 

Assume for a moment that the nucleon moves in the positive e^ direction, then p"*" 
and z~ are the only nonvanishing components of p and z, respectively. The infinite 
momentum frame can be visualised as the limit p'^ ^ Q ^ oo with fixed P ■ z = p- z ^ 1 
where Q is the large scale in the process. Expanding the matrix element in powers of 
l/p"*" introduces the power counting in Q. In this language, twist counts the suppression 
in powers of p"*". 

Similarly, the nucleon spinor N^{P^X) has to be decomposed in "large" and "small" 
components as 

N^{P, A) = (A+ + A-)A^^(P, A) = A^+(P, A) + N-{P, A) , (2.18) 
where we introduce two projection operators 

A+ = ^ , A- = ^ (2.19) 
2p-z 2p-z ^ ^ 

that project onto the "plus" and "minus" components of the spinor. Using the explicit 
expressions for N{P) it is easy to see that A+A^ = A^+ ~ Vp^ while A~A^ = A^~ ~ 
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2.3.2 Twist decomposition 

The twist decomposition of the nucleon-to-vacuum matrix element can be written in a 
form 

4(0|£"X(aix)4(a2x)rf;(a3x) |iV(P,A)) = F«({a,}, P • a;) W^'^V^'K (2.20) 

i 

where F^^^ G {5*1, 5*2, Pi, P2, Vi, . . . ,Vq, Ai, . . . , Aq, Ti, . . . , Tg} represent now nucleon dis- 
tribution amplitudes (DAs) and functions of definite twist: 



twist-3 
twist-4 
twist-5 
twist-6 



Vi, Ai, Ti 

T/„ A„ A„ T„ T„ TL 

(2.21) 



5i, Pi, V2, V^3, ^2, ^3, ^2, T3, T7 



'S'2, P2, V4, V5, A4, ^5, T4, T5, Tg 



The Dirac structures VF^^ and V^*"* can be read from Ref. [5H], and V^*"* contain the 
or projections of the nucleon spinor. 

The functions JF^*) and P'-*-' are related by 



Si = 


51 


2P 


X1S2 


= Si- 


^2 


Vi = 


Pi 


2P 


XV2 


= P2- 


Pi 


Vi = 




2P 


XV2 


= Vi- 




2V3 = 


= ^3 










Ai = 




2P 


XA2 


= -A, 


+ A2-A 


2^3 


= ^3 










ri = 


Ti 


2P 


XT2 


= Ti + 


T2 - 2T3 



(2.22) 



and each DA F^^ e Vj, Aj, Tj, Sj, Pj can be represented as 

P«({afc},P-x) = y"pxe-^^-^S.-'^"^P({x,}), (2.23) 

where the functions P*^*^ (xj) depend on the dimensionless variables Xj, < Xj < 1, = 
1 which correspond to the longitudinal momentum fractions carried by the quarks inside 
the nucleon. The integration measure is defined as 

= J dxidx2dx3 5{xi + X2 + X3 — 1) . (2.24) 

Analogously to (12.131) . the functions F^^\{xk}) poses the following symmetry properties 

r r)-f F^'K^2,x„x,) for P« G {V„T,} 
^ (Xi,X2,X3j - < ^ ^ A f _ ip(i) ^ rc /II ■ 12.25) 



as 



-P«(x2,Xi,X3) for P« G {5„P,-, A,} 
According to (I2.22p . the functions JF^*) can be written in terms of nucleon DAs P*^*^ 

.P«({aa, P ■ x) = ^-pl-^/(pW({a,}, P ■ x)) , n G {0, 1, 2} , (2.26) 



where / is a linear combination of F^^\ In addition, from fl2.11l) and Ref. [2H] one can 
read off the dependence of Dirac structures X^^^ and F*^*^ on x-coordinates: 

or or x^Xp . (2.27) 

Taking into account fl2.26l) and fl2.27l) the terms from (12.111) can be classified according to 



n = and X»F(*) ~ 1 
n = 1 and X(^)r(*) ~ X 
n = 2 and X^^^F^*) ~ x^x 



51, Pi, Vi, Vs, Ai, As, Ti, % 

52, V2, V2, V4, V5, A2, Ai, A^, T2, %, %, % (2.2^ 
Ve, A^, Tq, Tg 



After JF(*) functions are replaced by F^^^ and the Fourier transform (12.231) is employed, 
one ends up with corresponding three types of integrals which in LO take the form: 

2- /^'^ / / ^"^"^-^^^"^"^ ■ ■ ■] , (2.29) 

with kin being the momentum of the quark propagator. Similar but slightly more com- 
plicated integrals appear at NLO. 

The first integral in (I2.29P , corresponding to the first case in (12.281) , simplifies trivially 

to 



Vu (27r)V^)(g + - P n,a,)[F«({«fe}) 



(2vr)4 

J 

Pw[F«(K})---]M.=PE,-.a,-, , (2.30) 
It is convenient to introduce the notation 

A{ui,U2,Us)® B{ui,U2,Us) = I 'DuA{ui,U2,us)B{ui,U2,Us) . (2.31) 



and thus such contribution to the correlation function (12.51) can be expressed, both at 
LO and higher-orders, by a convolution 

= E-^?(^^}) ® ^^^HK}) • (2.32) 

i 

In order to evaluate the other two integrals from (I2.29P one employs 
y"rf^xa;«e'(''+'='"-P^i"^«^) " = ^(27r)^5(^)(g + fc;, - P ^ u,a,) j , 



Ujaj) j , (2.33) 
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as well as partial integratioij^, but we will leave out here further details. 



2.4 Derivation of LCSR 

Equating the correlation function (12. 5p results in several invariant functions that can 
be separated by the appropriate projections. Lorentz structures that are most useful 
for writing the LCSRs are usually those containing the maximum power of the large 
momentum p"*" ~ p ■ z. Hence, for the loffe current, we define the invariant functions A 
and B by 

z^A+T^ = {p-z){Mr,A+(^^l3}N+{P), (2.34) 
where A and B depend on the Lorentz-invariants = —q^ and P'^ = {P — q)'^. 

2.4.1 Correlation function versus form factors 

Next we relate the correlation function (12. 5p to nucleon form factors. The correlation 
function can be written as 

W <!) = \ p,2 E ^(0) ^)) (^(^'' ^)br(0)l^(^)> + ■ ■ ■ (2-35) 

^ s 

where the leading term is the nucleon contribution and the dots stand for higher reso- 
nances. 

Inserting (12. ip and the matrix element of the loffe current 

(0| r/(0) \N{P)) = Ai N{P) , (2.36) 

using Xls ^(-f =/' + ^N, and taking the projection suitable for our calculation, 

one obtains 

A+ = jji^^ {2(P ■ z) F,{Q') iV+(P) + {P ■ z) F,{Q') ^^N^P)}+--- 



(2.37) 



Employing partial integration in the second and third term of ()2.29p one gets 

-{ai-a„) / dui / dvi / dum 
^ Jo Ji Jo 

all possible l,m,n ^ {1, 2, 3} permutations for which a; — a,„ ^ 



and 



Vu e"*^^ {{uk})-, 



--(a;-a„) / dui / dvi / dwi / dum 

4 Jo Jl Jo Jo 

^g-^P■a.(«,(a,-a„)+n„(a„-a„)+a„)^(^)(|^^^ 1 - U;; - U„,}) , 

all possible l,m,n £ {1, 2, 3} permutations for which ai — a„i ^ 
Note that the surface terms that should vanish in Borel transformation have been already neglected. 
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Now, by comparing (12.371) and fl2.34p one gets 

-Fi(Q2) + ... = A{Q\P'- 

.F,{Q') + ... = E{Q\P''). (2.38) 



T7jn'^\^... - ^rn2 D/2^ 



2.4.2 Light-cone sum rules 

We can calculate A and B perturbatively in terms of quarks and gluons. 

Furthermore, (see, for example, Ref. [23]) we can formally write a dispersion relatioijf] 

1 r,l^B{Q',s) 



B{Q\P'') = - / ds / . (2.39) 





Now, by making use of the quark-hadron duality (analogously to, for example, Ref. 
[50] ) the effects of higher- resonances cancel between the left and right hand side of fl2.38p 
and one ends up with the sum rules 

2Ai ^,^2^ 1 r. ^raAiQ^s) 



1 

<Fi{Q^) = -J ds 



where sq is a convenient mass cut-off, which in our case corresponds to the continuum 
threshold taken at Roper resonance, Sq ~ (l.SGeV)^ [28], and which eliminates contri- 
butions other than nucleon (continuum subtraction). 

In practice, one can imagine to perform a power expansion of expression fl2.40p in 
the variable P'^ . To improve the convergence of this expansion one then employs Borel 
transformation 



^^^Ml\F{X)\= lim ^ ^ 



F{X) 



dX 



(2.41) 

\X\=nMl,X-^co 



Tin 

and of particular use in further calculation is a Borel transformation (see, for example, 

Bx^M%[l/ix - X)] = e-^'^^BlMl . (2.42) 



^Here we use a non-subtracted dispersion relation for which the condition 

A 
B 

has to be satisfied. 



Um^^oo R (Q^, s) = 
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Hence, taking a Borel transformation of the left and right hand-side of Eqs. (I2.40p 
in variables X = P'^, one ends up with the sum rules 

F2(Q2) = ^ dse^~'^^'"^^/^BimB{Q\s). (2.43) 
Aivr Jo 

For the Borel mass Mb which can be viewed as a matching scale of hadronic and partonic 
part of the calculation, we, as in Ref. [2H], take = 2 GeV^. 

We note here that in this draft we adopt the derivation of the sum rules based on 
the dispersion relations (12.391) and Borel transformation (I2.42p which then lead to (I2.43p . 
After calculating the needed imaginary parts and performing the necessary integrals it 
can be shown that this approach is equivalent to the one employed in Ref. [28j but 
perhaps more suitable for NLO calculations. 

The imaginary parts of the functions of interest are listed in App. [Dl 



2.5 Mtv dependence 

Due to calculational difficulties at NLO order, it is convenient to perform an expansion 
in mass Mn, i.e., schematically, the contributions to the correlation function can be 
presented by 

a + oM^ a a 

The first term corresponds to taking Mn = 0, while the first two terms correspond 
to taking = while retaining Mn proportional terms in the calculation. In calcu- 
lating NLO contribution we will try to investigate these two cases. In loop calculations 
additional In(M^) terms appear, and thus in the = approximation coUinear sin- 
gularities. 

The M^r-dependence of the general Lorentz decomposition of the nucleon matrix 
element (12.111) can be summarized as 



M^proportional X'-'^Y^^ 
M^-proportional X'-'^Y^'^ 
M^-proportional X^^^F^*) 
M|r-proportional X^^^F^*) 



Vi, Ai, Ti 

Si, Pi, V2, V3, A2, A3, T2, T3, % 

<S2, V2, V4, V5, A4,, A5, %, %, % 



(2.45) 



Obviously for M^v = only the terms proportional to JF^*) G {Vi, Ai,Ti} will contribute 
to T^. But for Mat 7^ the contributions proportional to J-'^'^^ e {Vi, Ai,Ti} will contain 
also M^-proportional (n 7^ 0) terms. 



3 LO contributions 

As a preparation and necessary ingredient of the NLO calculation, in this section we dis- 
cuss in detail LO contributions to correlator function (12. 5p . We present the complete 
twist-3 and twist-4 results and devote the last part of the section to the analysis and 
discussion of gauge invariance. 
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3.1 General structure and properties of LO contributions 

Using general Lorentz decomposition of the nucleon matrix element of the three-quark 
operator given in Eq. (12. lip and interpolating nucleon current of the form (12.71) . one 
obtains the general form of the LO contribution to the correlation function (12.51) 



T, 



LO 



.1 f f ^!^hlLpikinx\^ 

47 ^"^^ 7 (27r)4' ^ 



(27r)4 

X < e„ J^^'\P -x, 0,0) Tr 
+e„ J^(*)(0,P-x,0) Tr 
+ed J^«(0,0,P-x) Tr 



nit 



(3.1) 



and in following we neglect quark masses m„ = rrid = 0. Taking into account the form of 
it is easy to see that only the terms proportional to JF^*) g {Vj, contribute if loffe 
current, i.e., Fi = C'jx and r2 = 757'^, is taken for interpolating nucleon current. Namely, 
the matrices X« (see Ref. |28]) proportional to J^^*) G {Vj,Aj} and J^^'^ G {Sj,Vj,Tj} 
consist of odd and even number of 7 matrices, respectively, and hence the traces from 
(13. ip vanish in latter case. 

The contribution corresponding to the first case of (I2.28P takes the form 



T. 



LO 



X |e„ Tr 
+eu Tr 
+ed Tr 



(x«)^ri 



and using (I2.32p we can write 



(3.2) 



(3.3) 



where from comparison with (13.20 the definition of M^'^'^'^\{uk]) is obvious. As ex- 
plained in Sec. 12.3.21 for the other two cases listed in (I2.28P slightly more involved 
formulas appear. 

For loffe current (12. 8p . there are two twist-3 contributions (P^*^ G {Vi, ^i}) calculable 
by (13. 2p with X*^*) = 1, and steaming from 









Vi = 






(75iV). 


A = 









(3.4) 
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Figure 2: LO diagrams contributing to the correlation function 02.51] . On diagram B the 
nomenclature for NLO diagrams is sketched. 



Furthermore, when using loffe current, there are the two twist-4 contributions (F*^*) G 
{VJj,743}) calculable using fl3.2p with A^'^*) = 1/2, and two twist-4 contributions (F*^*) G 
{V2, A2}) corresponding to the more involved second case of fl2.28p . These contributions 
stem from 









V3 = 


1/2 V^3 










^3 = 


1/2 A3 










V2 = 


1/(2P. 


x){Vr~ 


V2 - V3) 






A2 = 


1/(2P 


x){-Ai 


+ A2- A3) 







Our main interest in this work are LO and NLO contributions to twist-3 contributions 
and the generalization to other first case contributions. Hence, we do not discuss in detail 
the calculation of second and third case LO contributions from fl2.28l) . nor higher-twist 
contributions, but rather refer to [2H] and references therein. 

In order to simplify higher-order calculation it is convenient to extract the Feynman 
rules and perform the calculation in momentum space. The three terms contributing to 
(13. ip correspond actually to three Feynman diagrams presented on Fig. [2l 

The Feynman rules for the contributions corresponding to first case of (12.280 are 
listed in Appendix |X1 Using these rules it is trivial to write down the contributions 
corresponding to the LO diagrams in Fig. O As it should, these contributions agree 
with corresponding terms in Eq. (13. 2p and 

= • (3-6) 
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Table 1: LO twist-3 contributions corresponding to the diagrams of Fig. [21 The contri- 
butions M.^'^^^ can be obtained using (13.71) . 



A y<A VI 


p 


K ACV^ 
^CAl 






Taking into account (12.121) it is easy to see that 



M 



n2,Mi,M3) for F«G{l^i,\/3} 

M2,til,^^3) for F«G{Ai,A3} 



(3.7) 



It is also obvious thal0 A^^'^i({Mfc}) = M^^^^Hu^}) = 0. 
Hence for the first case of (12.281) one can write 



-M^O'«(K}) 



4 

X je„Tr 



+e,Tr [X«(C7 



V(niP-g)2 {u2P-qY 



7^X«(C7a) 



Aj 757 



A 



(n3P - g)2 



where "+" sign in the first term corresponds to F^^^ E {Vi, V3} and "— " sign to F*^*^ G 
{Ai, A3}. Finally, note that, since we are interested in NLO calculation in which dimen- 
sional regularization will be employed, one has to use D = 4 — 2e dimensions also at 
LO. 



3.2 Twist-3 results 

In Tabled] we present the LO twist-3 contributions to diagrams displayed on Fig. [21 The 
complete LO contributions to amplitudes A4]f'^^^ are given by 

M]f'^'\ui, U2, Us) = M^'^'^im, U2, Us) + M^'^'\ui, U2, Us) + M''/'\U,,U2, Us) 
= M^'^-'^Ui, U2, Us) ± M^'^'\U2, Ml, Us) + Mf^'\ui, U2, Us) , 

(3.9) 

^ This result is expected since Ai and A3 are antisymmetric in ui ^ U2, while A^'^'(') (tVI^^)'-^'^')) 
contribution is obviously "blind", i.e., symmetric, to ui ^-s- U2 exchange since the photon couples to 
d-quark (which carries the momentum fraction U3 - see Fig. [5]). This property no longer holds at NLO 
order at which gluon can couple to u-quarks (see Fig. [3]). 
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where "+" sign in the second hne corresponds to Vi contributions and "— " sign to Ai 
contributions. The correlation function is given by (13.31) and taking into account the DA 
symmetry properties (12.251) . one obtains 

{M^'^'\ui, U2, U3) ± M^'^'\U2, Ui, M3)) ® F^'\uu U2, Us) 

= 2M^'^'\ui,U2,U3)®F^'\ui,U2,U3), (3.10) 

but it is advantageous for NLO calculation to leave full ui and U2 dependence of the 
amplitude A4. 

Our LO as well as NLO results for can be presented in terms of six invariant 
functions which multiply the Lorentz structures P^M^, ^, q^Mj^, ^, 7^, and 
7^ (pAf^. Thus it is advantageous for future calculations to present the results in a form 

Mf^'^^W}) = c^X(M) P^^N + cj«({n.}) p, 4 

+C^f{{u,}) 7, + C^^fuMnk}) 7. . 

(3.11) 

The corresponding LO twist-3 coefficients are given in Table [21 

Finally, we present the results for the invariant functions A and B defined in (I2.34p . 
Multiplying with z'^A^ the invariant functions A and B are projected. It is easy to 
see that these correspond to the invariant functions multiplying Pfj^M^ and ^, and 
thus 

A'--^-\Q\{P-qf) = C]^X(K}) ® ® 

(3.12) 

B'---'-\Q\{P-qf) = C'fJ\{u,]) ® V,{{u,]) + C]Pf\{u,]) ® A,{{u,]). 

(3.13) 

The results presented here are in agreement with the results from [28] . 
3.3 Twist-4 results 

The twist-4 contributions corresponding to V3 and A3 are obtained analogously to the 
twist-3 results discussed in the preceding subsection. The LO coefficients corresponding 
to (13. lip are given in Table El Note the dependence on e = (4 — D)/2. 

Furthermore, without going into much detail, we present also the contributions cor- 
responding to V2,A2 - see (13. 5p . We therefore introduce the shorthand notation 

^123 = V,~V2- Vs, 

A123 = -Ai + A2- A3. (3.14) 
Furthermore, as in [28], we use the definition 

dVi J dUmF{{vi,U.m,l - Vl - Urn}) . (3.15) 
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Table 2: LO coefficients (13.111) corresponding to Vi and Ai DAs. 



C 



hO,Vl 



c 



'Loyi 



^Loyi 



(jhoyi 



^Loyi 



L0,A1 
P,,Mn 



2uieu 





+ 2uiP ■ q 












+ 2uiP ■ q 


-ulMl 









+ {ui U2) 

+ {Ui ^ U2) ) + 



Q2 + 2miP ■ q - u\M 



2. .2 +{Ui^U2) ] + 



TV 



g2 + 2u^P ■ q - ulMl 
Q2 + ■ q - ulMl 







e^{P-q- UiM\ 



N) 



g2 + 2uiP ■ q - ulMl 



+ {ui ^ M2) ) + 



erf(-2P • q + 



g2 + 2u-iP ■ q - ulMl 



Q2 + 2usP ■ q - ulMl 



C 



^LO,Al 
^P^<i 



c. 



c. 



LO,Al 
LO,Al 



(^LO,Al 



^LO,Al 





+ 2uiP 


■ g- 






— 








+ 2uiP 


•g- 


ulMl 




-P 








+ 2uiP 


■ g- 






e 








+ 2uiP 


■q- 


u\Ml 



{Ui ^ U2) 



{Ui ^ U2) 
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Table 3: LO coefficients (13.111) corresponding to V3 and ^3 DAs. 



C 



LO,V3 



(3 - t)uieu ~ TV" -2(1 - e)M3erf 

+ (Ml ^ ^2) + 



+ 2uiP ■ q - ujM], 



g2 + 2U3P ■ q - ulMl 



^Loys 



^Loys 



(3 - e)eu 





+ 2uiP ■ q 


-ulM% 




UiMle 






+ 2uiP ■ q 


-ulMl 




+ 2miP ■ g 


-u\Ml 




Mi(l -e; 






+ 2uiP ■ q 


-u\Ml 









+ {Ui ^ M2) ) + 



2(l-e)erf 



Q2 + 2U3P ■ g - ulMl 



+ (Mi ^ U2)j + 
+ (Mi ^ U2) ) + 



g2 + 2M3P ■ q - ulMl 

ecd 

Q2 + 2M3P • q - ulMl 



c 



LO,A3 



^LO,A3 



Q2 + 2miP ■ g - ulMl 



C. 



LO,A3 







{ui ^ U2) 



- {Ui ^ U2) 



^LO,A3 



-(1 - e)uiMleu , s 
Q-^2u,P.q-u\Ml - ^ 



(1 - e)e„ 



g2 + 2miP • q - u\Ml 



(mi ^ M2) 
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Here F is the nucleoli DA or the combination of nucleon DAs that depends on three 
valence quark momentum fractions, and the integration over one momentum fraction 
has already been performed using 5(1 — vi — Um — Uk). Note that in this notation, which 
follows closely Ref. p8], F{ui) is not a simple function of ui and that the form of the 
function itself depends on ui = Ui, U2 or U3, i.e., whether ui corresponds to the momentum 
fraction of the first u-quark, second u-quark or d-quark^. Note that F(ui) = ±F{u2) for 
F{ui,U2,U3) = ±F(-U2, Ml, W3), respectively. 

The function F{ui) depends on only one momentum fraction and enters the expression 
analogous to fl2.32p . but that expression contains the integration with amplitude A^(Mj) 
only over one remaining momentum fraction Ui. This property holds in LO (where the 
dependence on mi, U2 and is clearly separated) while in NLO the expressions will 
be more involved. For tabulated LO contributions see Table HI The contributions to 
functions A and B are given by (1E.3P and (lE.Sp . 



3.4 Gauge invar iance 

Due to the presence of nucleon interpolation current, the condition of gauge invariance 
takes for the correlator function (12.51) the form 

q,T^ = Qn (Oh(0) |iV(P)) , (3.16) 

with Qr^ = 1 for the proton interpolation current and full electromagnetic current (12. 6p . 
while when one considers e^- and e^-proportional parts of the electromagnetic current 
separately Qn = 2e„ = 4/3 and Qn = = —1/3, respectively (in the neutron 
usual, ^ ed)- Furthermore for loffe current we have 

(Oh(0)|iV(P)) = AiM;vAr(P). (3.17) 

In the preceding subsections we have presented the complete LO twist-3 and twist-4 
contributions to the correlator function, i.e., in contrast to Ref. [2H], not just the terms 



^Hence, the shorthand expression (|3.15p encompasses three functions 

dvi J du2 F{vi,U2,l ~ vi - U2), 

dvi / du3 F{vi,l - vi - U3,U3), 
Jo 

dv2 J dui F{ui,V2,l - ui - V2), 

U2 rl — V2 

dv2 I du3F{l - V2 - U3,V2,U3), 



1 "'0 
^ ('tis rl — V3 

F{u3) = / dvs duiF{ui,l-ui-V3,V3) 
Jl Jo 

dV3 du2F{l~ U2- V3,U2,V3) . 

1 Jo 

Strictly speaking, it would be better that different names are introduced for these three functions instead 
of using argument to determine the form of the function. But for historical reasons and simplicity of 
the notation we adopt this notation hoping that it will not lead to too much confusion. 
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Table 4: LO coefficients analogous to (13.111) and corresponding to Vi23('Ui) and A^aiui 



LO,yi23 , 



{a 



hO,Vl2i I 



• e„((2 - e)Q^ + 2(1 - e)uiP ■ g - (1 - e)M>f 
(g2 + 2miP ■ q - ulM%f 

e,(-2e + 4(1 - 6)7x3^ ■ 9 + 2(e - 1)^^M^) 
(g2 + 2M3P . q - MiM2,)2 

' e„MiM^ 

. (g2 + 2uiP ■ q - ulM%Y 

eu{2P ■ q - UiM%) 
(g2 + 2miP ■ q - m2M2,)2 

■ -euMl 

. (g2 + 2u^P ■ q - ulMlY 

-CuMliQ^ + UiP ■ q) 
(g2 + 2^ZiP ■ q - ulMlY 



{ui ^ U2) 



(g2 + 2M3P • q - ulMlY 



. , ed(4P-g-2^i3M2,; 

(Ml ^ U2) 



(g2 + 2M3P • q - ulMiY 



{ui <^ U2) 



-2edM] 



N 



(g2 + 2n3P ■ q - ulM%Y 



{ui ^ U2) , 



N 



g2 + 2?i3P • q - ulMl 



eu{-P-q + UiM^ 



(g2 + 2?iiP ■ q - v\MlY 



{ui ^ M2) ,0 



• e„ (-eg2 + 2(1 - e)mP ■ g - (1 - e)n2M2,; 



(g2 + 2niP • q - ulM%Y 



, -{ui ^ M2) ,0 



-e„MiM 



AT 



^(g2 + 2niP-g-u2M2,)2 
eu{-2P ■ q + UiMl) 



(g2 + 2miP ■ q - ulM 



euM^ 



N 



-{ui ^ U2) ,0 



-{ui ^ M2) ,0 



-{ui ^ M2) ,0 



(g2 + 2miP ■ q - ulMlY 
-e^Ml{-tQ^ + (1 - 2e)^iP ■ g - (1 - e)^^M^; 

(g2 + 2miP ■ q - uiMiy 



'{ui ^ M2) ,0 



e,(P ■ g - uiM] 



N, 



(g2 + 2miP ■ q - ujM 



2A/f2 \2 ' 



-(Mi <^ M2) ,0 
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corresponding to the functions of interest A and B. This enables us to check the gauge 
invariance. By making use of the results given in Tables [21 El and H] one can easily see 
that for Mjf ^ the gauge invariance does not hold, i.e., that (13.161) is not satisfied, for 
separate ^*^*)-proportional terms nor for separate twists. 

Let us consider the expansion in given as outlined in Sec. 12.51 
In the M]\r = approximation, taking into account (13.171) . the gauge condition (I3.16P 
takes the simple form 

q>.T^ = 0. (3.18) 

In this approximation only twist-3 contributions proportional to Vi and Ai exist. As can 
be easily seen from LO results listed in Table Ej these contributions are separately gauge 
invariant, as well as, and parts separately. The same condition applies and holds 
for NLQj and higher-order contributions. 

Furthermore, we have checked the gauge invariance of the LO results in the = 
but Mn 7^ approximation to which Vi, V3, V2, Ai, A3, and A2 DAs contribute. This 
approximation corresponds to the first two terms in Eq. (I2.44p . Using the LO results 
given in Tables [21 [3] and [Hand taking = 0, one can show that 

q^T^ = Q^,XiMmN{P) (3.19) 

is satisfied when the sum of all contributing terms is taken into account (i.e., both twist-3 
and twist-4 contributions) and the asymptotic forms of twist-3 DA^ are used. There 
are no, at least at this order of expansion in M^v, conditions on twist-4 DAs. 

Hence, gauge invariance can be satisfied order by order in the expansion in (I2.44p 
with possibly some additional conditions on the form of DAs. For the check of higher 
order terms in Mn one should calculate the complete LO higher-twist contributions to 
the correlation function. 

4 NLO contributions 

We are finally ready to address in this section the NLO contributions to the correlator 
function (12.51) . We will present the results obtained in M^v = approximation (cor- 
responding to the first term in (I2.44p ). discuss the problems encountered in Mjv 7^ but 
Mff = approximation (corresponding to the first two terms in (I2.44p ). and outline the 
obstacles present in general Mjj 7^ calculation. 

4.1 Topological structure 

Three LO diagrams displayed in Fig. [2] lead to 30 NLO diagrams. The nomenclature 
we use can be deduced from Fig. [21 (diagram B): at NLO the gluon is attached in all 
possible ways. Typical NLO diagrams are presented in Fig. [3l 

The contributions of Bij diagrams one can obtain from the contribution of Aij dia- 
grams by Ml U2 exchange analogous to (13. 7p . The similar relation connects diagrams 
C24 and C23, as well as, C14 and C13. Taking all this into account, there are 8 more 

^ We will use this condition of gauge invariance to check the NLO results and resolve the 75-ambiguity. 
''if one, as natural, demands gauge invariance of e„ and terms separately then = 1/3 and 
A" = is enforced, while for the sum of e„ and terms = is sufficient. 
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complicated diagrams to calculate {A12, A23, A24:, A13, AU, C12, C24, C14) and the 
rest are either proportional to LO (Xll, X33, X44, X22, X34) or obtainable from above 
mentioned symmetries. 

Our calculation is performed in Feynman gauge. The colour factors were calculated 
using usual SU{Nc = 3)-algebra relations. In the first group of diagrams where gluon 
couples to the same quark line the colour factor is Cp = {N^ — l)/{2Nc) = 4/3 while 
in the second in which gluon connects two different lines colour factors equal {—Cb) = 
-{Nc + l)i2Nc) = -2/3. As mentioned in App. \M since we are describing nucleon as 
colour singlet state of three quarks (e"*'^), the Nc = 3 choice is enforced. 



4.2 Using dimensional regular izat ion and resolving 75 ambigu- 
ity 

In the following NLO calculation we take 

P'^ = Mjf = (4.1) 

and consider diagrams with massless on-shell external legs. Hence, apart from UV diver- 
gences the IR divergences of collinear type appear (there are no "true" IR divergences). 
Both divergences are regularized using dimensional regularization in 

D = A-2e 

dimensions. 

We introduce the abbreviations 

Tuvie) = ne f^^'^l^^^l^ i^^^T = 7 " 7 + ln(47r) + 0(e) , (4.2a) 

TiRie) = T{l + e)^^^^^^^{Anr = —^+^-HAn) + 0{e). (4.2b) 

The first F function on the right-hand side of Eqs. 04.21] originates from the loop 
momentum integration, while the integration over Feynman parameters produces Fs 
collected in a fraction. Consequently, the singularity contained in F(e) appearing in 
fl4.2al) is of UV origin, while the singularity contained in F(— e) appearing in (]4.2bp is of 



infrared (IR), i.e., collinear, origiiu. From F(2;)F(1 — z) = it/ sinvrz one can see that to 
all orders in e 

Tuvie) = -TjR{e). (4.3) 

Nevertheless, we find it useful to keep track of the origin of the UV and collinear singu- 
larities (for details, see also Refs. [52l [53]). 

In dimensional regularization the "trivial" self-energy diagrams (Xll, X33, X44), as 
well as, "trivial" 3-point integral (X34) vanish if one allows that UV and IR divergences 
to cancel. However, as mentioned above, we adopt here an approach of consistent track- 
ing of UV and collinear singularities, and their separate removal by renormalization and 
factorization, respectively. 



^The UV divergent integrals are finite in _D < 4 dimensions, while IR ones in Z) > 4 dimensions. 
Since we regularize both in D 4 — 2e dimensions, r(e) represents a signature of UV divergence, and 
r(-e) of the IR one. 
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All 



A34 




A13 A14 



Figure 3: Typical NLO diagrams. The contributions of "self-energy" diagrams A33 and 
^444 is equal to the contribution of diagram All. Diagrams Bij and Cij are similar to 
the ones presented above, while the nomenclature is sketched on Fig. [2l 
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Additionally, when calculating the contributions to corresponding to Ai distri- 
bution amplitudes we encounter 75-ambiguity - see App. [B] for details. In these cases 
the general Lorentz decomposition fl2.1ip and the choice of nucleon interpolating current 
(12. 7p lead to the appearance of the traces with one 75 matrix. At NLO these traces con- 
tain contracting 7 matrices and, as such, trace 75-ambiguity. Moreover, after the trace 
operation is performed one is left with one or more Levi-Civita tensors which get con- 
tracted with additional 7 matrices. Hence, the ambiguity related to the use of Chisholm 
identity is also present. 

We choose to use the naive-75 scheme |54j- We could choose to use HV scheme 
[55| [56] but then we would have to know or somehow calculate the terms which remove 
"spurious" anomalies violating Ward identities. Moreover we would have to use HV 
scheme also for the calculation of otherwise nonproblematic contributions corresponding 
to Vi. 

We remember that the choice of general decomposition (12. lip is not unique and that 
using Fierz transformations one could get the representation in which there is no trace 
and as such no 75 ambiguity (no trace ambiguity and no appearance of Levi-Civita 
tensor). Hence, the intermediate appearance of the problems with 75 are caused by our 
choice of the Lorentz decomposition of the nucleon matrix element and by the choice 
of the interpolating current. One can, for example, use the Lorentz decomposition of 
the form X^^Ya '' (see App. A in Ref. [25] for useful relations) which when used with 
(12. 7p does not lead to the appearance of the trace. But the price to pay when using this 
decomposition are much larger expressions, e proportional terms at LO even for Vi and 
Ai etc.. 

Nevertheless, that possibility lead us to the correct way to handle contractions using 
Chisholm identity: "follow" the fermion line (as usual, that means to go opposite the 
fermion line) and always perform the contraction of the Levi-Civita with the "last" 7 
matrix (with an open index) on the d line. The generalized recipe follows that one should 
write also the traces as a part of an expression obtained following the fermion lines - 
remember that the essence of trace ambiguity is loosing the cyclicity of the trace - see 
App. El 

When calculating NLO contributions in Mjv = approximation, we have used the 
gauge invariance (13.160 as a check and a help to resolve 75-ambiguity. After adopting this 
simple recipe - to write all parts of expression "following" the fermion lines: opposite 
d-line, along u-line (C ■ ■ ■ present - see App. opposite u-line, and to perform 
all evaluations obeying that order - we obtain the gauge invariant NLO results as one 
should. 

4.3 Twist-3 and M^v = 

Let us first investigate the approximation in which we neglect nucleon mass completely 
and take consistently throughout the calculation M^v = 0. As can be seen from (I2.45p . 
the general Lorentz decomposition of nucleon matrix element of three quark operator 
(12. lip has for M^r = only three terms: the ones proportional to Vi, Ai, and 7^. The 
tensor contributions vanish for our choice of interpolating nucleon current, and we are 
left with two contributions convoluted with Vi and Ai. 

^mIm«=o = MZ^{ui,U2,U:i) (g) Vi(mi,M2,«3) + ^^^^^(mi, ^2, ^s) ® v4i(ni, M2, M3) • (4.4) 
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As explained in, for example, Ref. [2S|, since Vi and Ai have different symmetry 
properties (see fl2.25l) ). they can be combined together to define a single independent 
twist-3 nucleon distribution amplitude: 

$3(0:1, X2, X3) = Vi{Xi,X2, X3) - Ai(xi, X2, X3) . (4.5) 

Now, by taking into account that the convolution of symmetric and antisymmetric func- 
tions gives 0, as well as, relation (14.51) and symmetry properties (I2.25p . one can write 
(USD as 

+ \ i^^^i'^l^ ^2, U3) - M^^{U2, Ui, M3)) O Ai{ui, U2, U3) 

= \ "2, U3) + MI^{U2, Ml, U3) - M'^^iUi, U2, U3) + Mf{Ui, U2, U3)) 

(8)$3(Mi,M2,M3) 

= A1*(mi,U2,U3) O <I>3(mi,U2, U3) . (4.6) 

From (14. 6 p the definition of A^* is obvious. 

4.3.1 Renormalization and factorization of collinear singularities: $3 

We start here with the explanation of renormalization procedure for (14. 6 p in which 
is given in terms of convolution of only two functions A^*(mi, ^2, M3) and $3(^1, U2, W3). 
We follow closely Ref. [52] . 

The amplitude 7V(*(mi, U2, U3) is of the general form 

M = Zr 

where 

A^LO _ ao + eai + 0{e^) , (4.7b) 
= {Tuv{e)[b^'' + ebr + 0{e')]+Tjn{e)[bi^ + ebi^ + 0{e')]}^^''^^ 



^curr 



An 



(4.7a) 



Q 



(4.7c) 



and Ai^'~' and Ai^^'^ are calculated from LO and NLO diagrams from Figs. [2] and [3l 
respectively. 

The bare coupling constant can be defined in terms of the running coupling con- 
stant a;s(/i|.) a^ 

= c^sif^l) ( 1 - ^^Po-] ( 41 ^ [^r^^(e)]-^ , (4.8) 



Ait 



and to the order we are calculating this essentially means that the bare coupling is 
replaced by the renormalized one and no singularities are removed as yet. 



In this as in the rest of the presentation we prefer to retain all terms in the expansion over e. 
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The expansion of the amphtude M. takes the form 



An 



+ 



where 



ao + eai + 0{e^) , 



2 \ ^ 



(4.9a) 

(4.9b) 
(4.9c) 



In our case the coefficients of the 1/e poles of UV origin, i.e., b^^ , are removed by 
renormahzation of the nucleon interpolating current (loffe current in this calculation). 
For that purpose .^curr has been introduced and it is of the form 



■^curr — 1 ~ Q;s(/^iJ,l 



L/curr 



(4.10) 



where /ifj is a scale at which the nucleon interpolation current is renormalized. The 
change of the scale of the renormahzation constant is given by 



(4.11) 



Hence, Ai takes the form 



4tt 



where 



(4.12a) 



(4.12b) 



and 



bo^ - O'c^^^j. ao 



ai 



For 



(4.12c) 
(4.13) 



the 1/e poles of UV origin vanish, and as only signature of their existence the logarithms 
Ccurr o-o ^^if^^R i/Q^) will remain in the end result. Notice that we have shown here that in 
principle coupling constant renormahzation and the renormahzation of the current can 
be performed at different scales, /^.^ and Generally, we write the amplitude A4 as 

an expansion in q;s(/x|j) and independent of /ij^. The truncation of this series in actual 
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calculation will introduce the dependence of the results on (frequently discussed in 
the literature). 

The remaining collinear singularities get cancelled by the renormalization constant 
of the nucleon distribution amplitude. Namely, 



= Mf^iui, U2, U-s) ® Z^{ui, U2, Us; Xi, X2, X3; fip) <S) ^3{xi,X2, X3; ftp) 



Mf^{xi, X2, X3; fip) (g) ^3{xi, X2, X3; . 



The renormalization constant Z$ is of the form 



Z^{{uk}, {xk}; Hf) = 5{ui - xi)5{u2 - X2) - ^^(/ij 



(4.14) 



O(a^), (4.15) 



where is a leading term of the kernel of the evolution equation for twist-3 DA: 



d 

/i^^$3({Mfe};/^^) = V$(K},{x4;/i^) ® <l>3({xfc};/x^) 



and 



47r 



(4.16) 



(4.17) 



The kernel V^^ was given in Ref. [7] and confirmed, for example, by the calculation of 
anomalous dimensions in Ref. [57]. Here we present it in a convenient form 

(K},K}) 
3 



^Cf {^{U3-X3) 
+2Cb \ 5{U3 - X3 



■d{Xx - Ui) H d[X2 - U2) 



X\ X\ — U\ 
Ml 



"2, 



X2 X2 — U2 
6 



3 ^ 1 
3^2 



9(xi - Ml) H 6(x2 - U2) 

Xi X2 



where hi = —hi is a helicity of quark i, while 



hlh2 ,(3^1 

^'3^2 



Ui + U2 



, (4.18) 



{F(K}, {xk})}^ = F{{uk}, {xk}) - 6{u, - xi)6{u2 - X2) J VzF{{zk}, {xk}) . (4.19) 
We stress that one should take Cf = '^Cb = 4/3. Furthermore, in fl4.18p one has to take 



i.e., the helicities of the quarks correspond to 

«Ta)«i(T)^Ta) ' 

which is in agreement with App.C from Ref. [28] . 
For the finite amplitude M. one then gets 



hz 







(4.20) 



(4.21) 



(NLO/ 



47r 



(4.22a) 
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where 



(4.22b) 



1 



1 

+ 

— e 



2\ ^ 



X 



2 \ f 



(4.22c) 



where /ii? is a (usual) factorization scale and the collinear singularities cancel for 

- ao ® V^^^ = . (4.23) 
With all singularities cancelled, we can now take the limit e — > and finally obtain 



47T 



-M''^^(K};/^^,/i^j,i) + 



with 



A^^°({xa) 



ao 



(4.24a) 



(4.24b) 



A<^^°({x.}; /.I, /.?,,) = {br - C^^Lai) - - a, ® F^^)) 
For loffe current 

nil) _ _ 2 

Furthermore, one can check in Table [1] that for our twist-3 results there are no e propor- 
tional LO contribution (in contrast to, for example, twist-4 contributions) and oi = 0. 
Taking this into account along with the conditions for cancelling UV 04.131) and collinear 
singularities fl4.23p . simplifies the NLO result fl4.24cp to 



(4.24c) 



(4.25) 



M^'^iix,}; /i^, f^i,) = br - b{^ + Hf^iim - c H^^Fm ■ (4.26) 

4.3.2 Renormalization and factorization of collinear singularities: Vi and Ai 



Although (14. 6 p and the analysis given in the preceding subsection are sufficient for ob- 
taining twist-3 results, we now turn to renormalization procedure for expressed in 
terms of Vi and Ai. i.e, as in (14.41) : 

^/^lAfjv=0 = -^^^("1' "2, Us) (g) Vi{ui, U2, Us) + M^^{ui, U2, U3) ® Ai{ui, U2, U3) . 

The crucial difference in comparison to the preceding subsection is that is no longer 
expressed as just one convolution but rather as a sum of convolutions. As we shall see. 
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there is a mixing between these terms. The procedure is similar to the one used in Ref. 
[SH] and, although the finite results are the same as the ones in preceding subsection, the 
experience that we gain in this subsection should be very useful for the Mj^ ^ case. 
It is instructive to write (14. 4p in a matrix form: 

= {Ml\u,, «3), -Mf u,)) ® ( ) • (4.27) 

Both and M^^ can be expanded as in (14. 7p and the UV renormalization of these 

expressions proceeds the same way as explained in the previous subsection. One ends 
up with the UV-finite expressions of the form (I4.12p : 



Mt\{u,}) = Ml^'^\{u,}) + ^^M^'^^'^^iu,}) + . . . , (4.28) 



and the conditions (14.131) have to be satisfied, i.e.. 



,UV,V1 _ ^(1) VI _ Q 

br'^^'-Cila^' = 0. (4.29) 

In order to cancel remaining collinear singularities one has to know the evolution 
kernels, i.e., renormalization constants, for Vi and Ai distributions amplitudes. One can 
probably derive it by additional one-loop calculation, or, as we will do here, make use 
of our knowledge of V$. Knowing the sjTumetry properties of Vi and Ai distribution 
amplitudes, we write Vi as 

V<s,{{uk}; {xk}; fi"^) = Vvi,vi{{uk}; {xk}; fi^) + Vvi,Ai{{uk}; {xk}; fi'^) 

+ VAi,vii{uk}] {xk}; A^^) + VAi,Ai{{uk}; {xk}] Ai^) , (4.30a) 

where 

Vvi,vi{{uk}] {xk};n^) = ^ {V<s,{ui,U2,U3]Xi,X2,X3; fi^) + V^{ui,U2,U3;x2,Xi,X3; fi^) 

+ V^{U2, Ml, Ms; Xi, X2, X3] fX^) + V^{u2, Ui^U^] X2, Xi, X^] /i^)) , 
yVl,Al{{Uk}; {Xk};fJ'^) = ^ {V^{Ui,U2,U3;Xi,X2,X3; f?) - V^{Ui,U2,U3;X2,Xi,X3; fl^) 

+ V$(U2, Ml, Ms; Xi, X2, X3; /i^) - 14(^2, Mi, Ms; X2, Xi, X3; /i^)) , 

VAi,vi{{uk}; {xfc};/i^) = i (V<i,(Mi,M2,M3;a;i,X2,xs;Ai^) + Vi(Mi, M2, M3; X2, xi, xs; /i^) 

-14(m2, Ml, Ms; Xi, X2, X3; /i^) - Vi(M2, Mi, M3; X2, Xi, Xs; /i^)) , 

VAi,Ai{{uk}; {xfe};/i^) = i (Vci,(Mi,M2,M3;xi,a;2,xs;/i^) - V$(mi, M2, M3; 0:2, xi, X3; /i^) 

-Kj,(m2, Ml, Ms; Xi, X2, X3; IX^) + Vi(M2, Mi, Ms; X2, Xi, xs; /i^)) . 

(4.30b) 
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Obviously, 

yvi,vi{{uk}', {xk}', jJ^) I symmetric in ui ^ U2 and symmetric in xi ^ X2 

H'i,Ai({wfc}; {xfc}; yU^) I symmetric in ui ^ U2 and antisymmetric in xi ^ X2 

VAi,vi{{uk}', {xk}', At^) I antisymmetric in Ui U2 and symmetric in Xi ^ X2 

Vai Ai{{uk}', {xk}', /i^) I antisymmetric in ui ^ U2 and antisymmetric in xi ^ X2 

(4.31) 

We can now substitute fl4.5l) and fl4.30al) in the evolution equation (14.161) and taking 
into account the symmetry properties with respect to Xk one gets 

d 

f^^-gj^ [^1 - ^1] {{uk}; /i^) = [Vviyi + VVi.Ai + Vai,vi + Vai,ai] {{Uk}] {Xk}] /i^) 

®[Vi-Ai] {{xk};fi^) 

= [Vviyi + Vai,vi] {{uk}; {xk}; /i^) ® Vi{{xk}; /i^) 

- [Vvi,ai + Vai,ai\ {{uk}; {xk}] /i^) O Ai{{xk}; /i^) • 

(4.32) 

Furthermore, the symmetry properties with respect to Uk allows us to write the evolution 
equation in a matrix form as 



(4.33) 



The DAs Vi and Ai obviously mix under renormalization and we can write 

'] (K}) = Z({n,}; {a;,}; fi') (^^ ] ({a;,}; /i^) , (4.34) 



where 



and 



z(K}; {xkh /i') = 1 - ^7^— v«(K}; {^4) , (4.35) 

47r — e 



T/(l) _T/{1) 
'l/l 1/1 ^1/ 



v(^H{«4; {^J) = Z'S JiV'''' ({^^5 {^'^i) • (4-36) 

\ ^Al,Vl ^A1,A1 / 

By substituting (gJl in (|42ZD, we get 
TJm,=o = (-M^'-^f ) (Kl) ® Z(K},{a;,};/x2) ® ( J (K};/x') 

= {M]:\Mf) {{xkhfi') «^^>'/^')- (4.37) 

The condition for cancelling coUinear singularities (14.231) now takes the more involved 
form 

C'^' + aj;' ® V^;;^i - «o ' ® = 0. (4.38) 
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Finally, we take e 



limit and obtain 




47r 



47r 



(4.39) 



with 




a, 



a, 



VI 
' 



Al 




(4.40) 



In our calculation = and taking this into account along with the conditions 

of cancellation of UV fl4.29p and coUinear singularities (14.381) . one finally gets 



and when one uses the actual values for biS one can see the agreement with (I4.26p . 
4.3.3 Results 

The cancellation of singularities for the M^v = case has been checked and shown for 
two equivalent representations: one corresponding to $3 and the other corresponding to 
Vi and Al DAs. In the latter case the mixing appears. 

In App. Owe list our finite NLO results contributing to the function of interest B. 
The function A cannot be accessed in = approximation - see (12.341) . For complete- 
ness sake we list both the results corresponding to Vi and Ai distribution amplitudes, 
as well as, to $3. The latter results are shorter and actually used in further numerical 
calculation. 

4.4 Away from = approximation 

Next we consider the approximation 7^ while = 0, i.e., the approximation 
corresponding to the first two terms of the expansion (I2.44p . 

We have calculated the NLO contributions corresponding to Vi and Ai. The UV 
singularities get cancelled as in the previous subsection. In contrast, the coUinear singu- 
larities cannot be cancelled considering only the contributions corresponding to Vi and 
Al. For example, both Cp^'^^ and Cp^^]^^ parts proportional to are 0, while NLO 
counterparts are different from and contain coUinear singularities. Obviously, mixing 
between Vi and Ai alone cannot cancel these terms even if we knew the new kernels 
H^i,vi,^Vi,Ai etc. corresponding to the 7^ case. 



(4.41) 
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Hence, it follows that the mixing with twist-4 DAs V3, A3, and even maybe V2 and 
A2, should be taken into account in this = but M^v 7^ approximation. This is 
similar to the observation from Sec. 13.41 that the gauge invariant results are obtained 
not twist- by-twist but order- by-order in M^. The problem which we then encounter is 
that we do not know these corresponding new kernels and they play a role not only in 
cancelling the singularities but also change the finite parts (oi, i.e., e proportional LO 
parts, are not for V3 and ^3 terms). 

For example, for mixing just between Vi, Ai, V3 and ^3 we encounter the unknown 
kernel of the form 

/ Vviyi Vvi,Al Vvi,V3 H^l,A3 \ 

V^iyi Kai.ai VAiys Vai,A3 

yV3,Vl Vv/3,A1 yV3,V3 H'S.AS 
\ ^ASyi Ka3,A1 yA3,V3 Va3,A3 / 

By substituting the conditions for cancellation of UV and coUinear singularities, for 
the finite NLO contributions one gets 



(4.42) 



lUvyi 
^0 



^0 



^V3yi ^ "1 



A3 



(1) 



A3,yi 



iUV,Al 7 /-R,A1 I ,UV,A1 1 / 2 /^2\ 



IR,A1 
'0 



uvy3 



6{«'^^ + 67'-ln(/.^,/g^ 



tuvys 



^1/3! V3 + ^l^ ® ^A3,y3 



,A3 



^curr 1 



V3 



iUV,A3 7lB.,A3 , ^UV,A3^ / 2 /^2\ 7,/R,A3 , / 2 



+ ( «r ® ^y3,A3 + «f ® VaU) - C^rr 4 



A3 



(4.43) 



But we do not know Vy^^ and V\^^^ and hence we do not know how to calculate the 
finite terms 



(4.44) 



4.4.1 Open problems in higher-twist calculations 

If we take 7^ 0, there are no coUinear divergences whose cancellation we should 
take care of. But there are additional open problems one should solve before attacking 
higher-order higher-twist calculations. 

For example, there is an open problem in the calculation of the NLO contributions 
to second and third case defined in (12.281) . Let us for a moment go back to coordinate 
space. When calculating NLO contributions one encounters the matrix elements of the 
form 

{Q\e^'%l{x)u'^{z2)d'^^{Q)\N{P,X)) , (4.45) 
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and 

(0| e'^'^ulix)u'^{z,)d;{z^) |iV(P, A)) , (4.46) 

where Zi and Z2 come from the gluon coordinates. In contrast at LO only matrix elements 
of the form 

{0\e^"^uUx)4iO)d;iO)\NiP,X)) (4.47) 

appear. For first case of (12.281) these additional coordinates, not apparently proportional 
either to x-coordinate (photon coordinate) nor 0, do not pose a problem but in the 
second (and third) case it seems not to be clear how to identify the coefficients and x 
in (EIH). 

These and similar problems are postponed for future investigations. 



5 Light-cone sum rules 
5.1 LCSR for Mn = case 

In Mat = approximation we cannot asses the function A and thus the form factor Fi. 
The form factor F2{Q'^) is given in terms of B{Q'^, P'^) by the sum rule (12.431) 

F2{Q') = ^ dse-^/''^lmBM^=oiQ\s). (5.1) 

The function B calculated for Mjv = case and to NLO is given in App. O Note 
that the convenient dimensionless quantity 

(5.2, 

has been introduced and that the function B has been expressed through this new vari- 
able, B{Q'^,P''^) —>■ B'{Q'^,W). For simplicity sake, the same name for the function B 
has been retained. Thus in Eq. (15.11) the change of variables s ^ w = {s + Q'^)/Q'^ has 
to be performed. 

Furthermore, note that i?7-terms (?7 > and rj <^) originating from the Feynman 
rules for quark and gluon propagators were explicitly kept in resulting logarithm terms 
throughout perturbative calculation. The analogous terms in denominators can be easily 
recovered resulting in W W + ir]. Hence, as it turns out, in our calculation the 
sign infront W and irj is the same both in denominators and logarithm termJ^ . The 
imaginary part can now be determined using the expressions listed in App. |Dl 

Let us see this in more detail. As can be seen from App. [Cl the function Bmm=o can 
be expressed in a form of convolution 

BM^=oiQ\W) = ^J2'^B,Fi^\M^=oiM,W;fil/Q')®F^'\{xk};fil) (5.3) 



"'^'^In practice, one often suppresses irj terms during calculation and recovers them when the analytical 
continuation is needed. This approach can in some cases (when more complicated functions appear) be 
non-trivial and can even lead to mistakes. Actually, it is much simpler just to keep track of irj terms 
from Feynman diagrams to resulting higher order expressions. In our work we adopt that approach. 
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with F*^*) denoting, as before, nucleon DAs (F*^*) G {Vi, Ai}, or F^^^ = $3 when the sum 
has only one term - see Eqs. (IC.ip and (lC.2p ). 

The imaginary parts of p(i) ]^^^Q{{xk}, W; fJ.p/Q'^) determine the imaginary parts of 
i3Mjv=o(<5^5 W). Furthermore, the hard-scattering part Tig p(i) j^^^(^{{xk} , W; ii^/Q"^) can 
be conveniently expressed as a sum of the terms of general form 

9{{xk}. W) /({x,}; /i|/Q^ lilJQ\ lil) , (5.4) 

where only (yf-functions possess poles leading to imaginary parts. The imaginary part is 
then determined from the imaginary part of the function g{{xk}iW) and, when from 
four integrations present in (15.11) two are performed, one obtains the following rule for 
separate terms contributing to Tq and leading to separate terms contributing to F2: 

(so+Q")/Q" , , 

■ -2 /n^^ 



AlTT J Ji 



X lmg{{xk},w) 

X fiixk}; liMQ\ AJQ\ A) ^^'^{^4; /^f) 



1 



1 /.l-x 



dx, \ dxj e-(i-^»)« /(^''^^^^ 



Ai 

X g{{xi,Xj,l - Xi - Xj},Q'^,so,Ml) 

X fi{xi, Xj, 1 — Xi — Xj}; jip/ Q , Q , /ijj) 

X F^^({Xj, Xj,l Xi Xj^] flp) . 



(5.5) 



The selected g functions (see Eq. ( IC.SI) ) that appear in our LO and NLO calculation of B 
are listed in TableOalong with corresponding 5? functions which contribute to F2 as shown 
in (15.51) . This table along with Eq. (15.51) thus provides us with necessary substitution 
rules for calculation of F2 from perturbatively calculated results for B summarized in 
App. O 

The resulting nucleon form factor F2 calculated to NLO takes the form 

F2(Q2, M|; fil, fil ,) = F^'^iQ', Sq, ; ;x|) + F^'^'^iQ^, Sq, M|; fi% /i^ J , 

(5.6) 

where 



71 I 

+F2 [Q ,Sq,Mj^-^p) 
+F2 [Q ,So,Mj^;idp) In — 



(5.7) 



where, as defined before, is the virtuality of the photon probe, is the coupling 
constant renormalization scale, /i^ ^ is the renormalization scale at which the renormal- 
ization of the nucleon interpolation current has been performed (often taken the same 
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as but in principle an independent scale), fijp is the factorization scale at which the 
coUinear singularities corresponding to the nucleon DA factorize, Sq is the scale corre- 
sponding to the continuum subtraction in LCSRs, and Mg is a Borel mass, which can 
be regarded as a matching scale of hadronic and partonic part of the calculation. 



5.2 Nucleon DAs 



rs = In, 

03 = y/^^i- 



(5.10) 



We refer to App. B of Ref. [2H] for detailed account of nucleon distribution amplitudes 
and list here only the selected expressions. 

For Mn = case only the twist-3 DAs are relevant: 

Vii{xk},fx^) = 120x1X2X3 [0;i(/x') + 0^ (/i')(l- 3X3)] , 
Ai({Xfc},Ai^) = 120XiX2X3(x2 - Xi)0j(/i^) , (5.8) 

or, equivalentlj0 

= 12OXiX2X3[0°(/i')+0+(/i2)(l-3X3)-0j(/i')(x2-Xi)] . (5.9) 

Here 



The normalization obtained using QCD sum rules amounts to [28] 

/7v(/io = IGeV^) = (5.0 ± 0.5) ■ lO^^GeV^ . (5.11) 

Actually the normalization of twist-3, -4 and -5 DAs is determined by three dimensionful 
parameters /jv, Ai and A2 that are well known from the QCD sum rule literature and 
correspond to nucleon couplings to the existing three different three-quark local opera- 
tors. The numerical values of the other two normalization constants, obtained by QCD 
sum rules [2^, are 

Ai(/i^ = IGeV^) = -(2.7±0.9) ■ 10"^GeV^ (5.12) 

and A2(/io = IGeV^) = (5.4 ± 1.9) ■ 10~^GeV^. For the evolution of these parameters we 
refer to, for example, |25j. 

The shape of the twist-3 DAs, i.e., the deviation from the asymptotic form, is de- 
termined by dimensionless parameters Vf and A^, while three more parameters (/f , 
and ) appear in twist-4 and at higher twists. The values of these parameters and thus 

As explained in, for example, Ref. [25], Vi and Ai have different symmetry properties and they 
can be combined together to define a single independent twist-3 nucleon distribution amplitude $3. 
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Table 5: Table of the substitution rules corresponding to Eq. fl5.5l) for (?„ functions 



(n = 0, ■ ■ ■ , 12) defined in Eg. ^M)- 



gn{Xi,Xj, W) 


gnixi,Xj,Q'^,so,M'^) 


1 


1 


{XiW -1 + 17]) 


Xi 



ln(l - XjW - IT]) _^ / sq + Q^ \ 

1 , 1 r /(i-t)g2\ ■ 



ln2(l - x,W - irj) 1 
(1 — XiW — irj) Xi 



XiJi [ V XiMl J 



ln(l — XiW — irj) 
{W + irj) 


-i: 




ln(l — XiW — irj) 
{W + ir])'^ 






ln^(l - XiW - iri) 
{W + iri) 






\t?{1 - XiW -iri) 
{W + ir])'^ 


-2x, 
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ln(l — {xi + Xj)W — iri) 
(1 - XiW - if]) 



1 

Xi 



In ( a; 



1 r^^^ 



+— 

Xi 



dt 



+ In 



Xi 



(i-t) 



exp 



(l-t)Q- 



- 1 



ln^(l - {xi + Xj)W - irf) 
(1 — XiW — irj) 



1 

X.: 



In ( X 



-1 1 +ln' 



Xi 



TT 



-y-7r2(l-5(x,-)) + 2Li2 



2 /--^^ ln(t(l 



El) - I 



(Xj + Xj ) ( 1 Xj "q? ) 

{l-t)Q' 



Xi ./^ 



exp 



XiMl 



- 1 



ln(l — (xi + Xj)W — irj) 
{W + irj) 



at - exp ' 



t 



XiMl 



ln(l — (xj + Xj)W — if)) 



Xi 



at — exp 



t^ 



XiMl 



ln^(l - {xi + Xj)W - ir)) 
(W + irj) 



-2 



dt 



In {t{l + I 



- 1 



exp 



(1 - w 

XiMl 



ln^(l - (xj + Xj)W - ir]) 
(W + ir])^ 



-2x. 



.£0+2: 



In (t(l + |) - 1^ /(l-t)Q2 
at — ^ 77, exp ' 



^2 



Note: As expected, the CQjSG Xj = corresponds to the results from previous page. 
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the shape of DAs are controversiaL The older prediction from QCD sum rules, which is 
sometimes referred to as the Chernyak-Zhitnitsky-like (CZ-like) model [SHI amounts to 

V^^ = 0.23 ±0.03, 

= 0.38 ±0.15, (5.13) 

while = 0.40±0.05, = 0.07±0.05 and /| = 0.22±0.05. In Ref. ^8j (Eq. (42,[28j)) 
the following values (we refer to them as BLW parameters) were introduced 

Vf = 0.30, 

= 0.13, (5.14) 

while ff = 0.33, = 0.09 and = 0.25. Finally, for the asymptotic DAs 

V,' = 1/3, 

= 0, (5.15) 

and ff = 3/10, f]" = 1/10 and /| = 4/15. Thus, for asymptotic twist-3 DAs 

4 = In, 0+ = 03- = O, (5.16) 

leading to 

<^3{{Xk}, /i') = ViiiXk}, /i') = 120 Xi X2 X3 /iv(/i') . (5.17) 

6 Numerical results 

In order to obtain numerical predictions for the proton form factor F2 calculated at 
twist-3 to NLO order we use the results (IC.18tlC.21l) . twist-3 DA $3 defined in 

Sec. 15. 2[ recipe (15. 5p and Table The additional necessary numerical values listed in 
preceding subsections and taken from Ref. p8], are 

So = (1.5GeV)^ 
Ml = 2GeV^ (6.1) 

and 

Xiifxl = IGeV^) = -2.7 ■ 10^^ Q^y2 _ ^q ^) 

We will check the sensitivity of the results on the variation of M^. Alternatively, one 
can use the modified value of Ai calculated using the expression (B10,|25j) from Ref. 
where is corrected by a factor (1 ± 71/12a;s/7r) from Eq. (24) in Ref. 



^corr.^^ _ loGeV^) = ((-3.4) - (-3.2)) ■ lO'^ GeV^ . (6.3) 
The one-loop expression for reads 

"^^^'^ = /5oln/.VA^' ^'-'^ 

where jSo = 11 — 2/3n/, with nf = 3 being the number of flavours. For A we take the 
value A = 0.2 GeV. In this work we do not take into account the DA evolution i.e. we 
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neglect the evolution (see, for example, of fNi.l^'^) and of Ai(/i^). One takes that 
these effects are small, especially since in the most part of the numerical analysis we fix 
the relevant scale to 2 GeV^. 

The LO twist-3 and twist-4 results calculated for M^r 7^ are presented in App. [El 
while the higher-twist results and corrections we take from Ref. [25 • The parameters 
for the DAs are given in Sec. 15.21 while for the exact form of the higher-twist DAs we 
again refer to Ref. [2H]- 

It is convenient to normalize the results for F2 and Gm to dipole form factor, i.e., 
GM/if^pGn) where 

(l + g2/o.7iGeV2)2- ^^'^^ 
6.1 Mn = case: F2 at twist-3 to NLO order 

In Fig. UK we present the relative size of the NLO correction (15.71) taken at yU/; = yul' = 
= Q"^ compared to LO prediction, i.e., the ratio F2^'^ / in dependence on Q"^ 
and calculated for the asymptotic twist-3 DA (I5.17p . Note that for that choice of scales 
only the p^^^'^'^ part from (15.71) contributes. In Fig. |4)d the complete NLO prediction 
F2 (15. 6p normalized to fipGo is displayed. In Fig. IHwe test the sensitivity of the results 
on the choice of by displaying the results for the default choice (16. ip . as well as, for 
M| = 1 GeV^. 

From Fig. one can see that, for the asymptotic DA and = (^"p = fJ'jn = Q"^ the 
ratio F2^^ / F2^ , and thus the NLO correction in comparison to the LO prediction is 
quite large, being 60-70% for 1 < < 10 GeV^, and it increases with Q"^. By decreasing 
to 1 GeV^ the ratio drops at higher Q"^ only slightly, i.e., to 65%. 

The sensitivity of the LO (dot-dashed and dashed hues) and NLO (solid and long- 
dashed lines) predictions for F2 on Mg is illustrated in Fig. Hb. One can see that this 
effect is large (comparable to the change from LO to NLO). 

In the following we adopt the default choice (16. ip employed in Ref. [28]. 

In Fig. the change of scales and ^ is investigated. As can be seen by 

comparing Figs. and[5H, by taking fJ^p = fJ^ji i = M'^ and retaining = Q^, the ratio 
of NLO correction to the LO prediction is lowered to 69% - 44% for 1 < < 10 GeV^ 
and it decreases with Q^. Thus, one can see that the ^^^^^'^^ and jr^^'^^^^ terms from 
Eq. (15.70 decrease the size of NLO correction. By changing /i^ to the ratio gets 
bigger since then there is no suppression due to the running of the a^. From Fig. 03 
one can see that the change of scales does not influence much the value of the complete 
F2 NLO prediction. In further calculation, if not specified differently, we, following [33] . 
take n% = nj, = = Ml. 

Finally, in Fig. [6] we investigate the size of the NLO correction and of the complete 
NLO prediction in dependence on the choice of DA (15.90 . We employ the asymptotic 
(EISD DA (solid line), BLW fimp DA (dash-dotted line), and the CZ-hke DA fCT|) 
(dashed line). In Fig. [6b the LO predictions are denoted by thin lines and NLO predic- 
tions by thick lines. From Fig. [6|l one can see that for 1 < < 10 GeV^ the ratio of 
NLO corrections to LO predictions amounts to 57% - 62% for the asymptotic DA, 61% 
- 76% for the BLW DA, and 66% - 86% for the CZ-like DA. Both LO predictions and 
NLO corrections are larger for the two DAs whose form differ from the asymptotic one. 

In conclusion, we state that the NLO corrections to F2 calculated at twist-3 taking 
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asymptotic DA, iur^ =iJ-f^ =I^r,i^ asymptotic DA, Hr"^ -iip^ -Hr^i^ =(f 




02468 10 02468 10 



[GeV^] [GeV^] 

(a) (b) 

Figure 4: NLO calculation of twist-3 contribution to the proton form factor F2 obtained 
for Mtv = (15 .(ins .71) . using asymptotic DA flHTTD and = n], = = Q^. ^he 

ratio F^^^'^^/F^*-* (i.e., the ratio of NLO correction to the LO prediction) evaluated for 
the defauh choice (El]) M| = 2 GeV^ (solid line), and for the test choice M| = 1 GeV^ 
(dashed line), b) The twist-3 prediction to F2 normalized to ^pGr,: NLO (thick solid 
line) and LO (thick dot-dashed line) for M| = 2 GeV^, NLO (thin long dashed line) and 
LO (thin short dashed line) for = 1 GeV^. 




02468 10 02468 10 



[GeV=^] [GeV=^] 

(a) (b) 

Figure 5: NLO calculation of twist-3 contribution to the proton form factor F2 obtained 
for = fl5.(iti5.7D and using asymptotic DA fl5Tfl) : a) The ratio Ff^^/F];^ (i.e., the 
ratio of NLO correction calculated at to the LO prediction) calculated using = fij^ = 
M^, while = (solid line), and = (dashed line), b) The twist-3 prediction 
to F2 normalized to fipGf): NLO for = and ^% = fJ'Ri = (thin long dashed 
line), NLO for = and fij, = fi% -^ = M% (thick sohd line), NLO for ^\ = M| and 
fJ^'p = fJ'R I = M"^ (thin short dashed line), and LO (thick dot-dashed line). 
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Figure 6: NLO calculation of twist-3 contribution to the proton form factor F2 obtained 
for Mtv = (15.6115.71) using /^/j = = /i/j 1 = The displayed results correspond 

to asymptotic DAs (solid line), BLW DAs (dash-dotted line), and CZ-like DAs (dashed 
line), a) The ratio F2^^°/F2^°, i.e., the ratio of the NLO correction to the LO prediction, 
b) The twist-3 prediction to F2 normalized to fipGn. Thick lines: NLO prediction. Thin 
lines: LO prediction. 



Mn = are large, with F^^'^'-'/Fl"^ amounting to cca 60%, but varying for different DAs 
and depending on the choice of renormalization and factorization scales. The sensitivity 
of both LO and NLO corrections on the choice of Mb is large and in the following 
we take the value (16. ip from Ref [2H]- In contrast to the dependence of F^^'^/F2^°, 
the dependence of the complete NLO prediction of F2 on the choice of renormalization 
and factorization scales is small and, if not stated otherwise, we take in the following 
fiji = fi'^p = fi%i = Mq. The results for different DAs differ significantly. 



6.2 Various contributions to Fi and F2 

In preceding section we have discussed the size of the NLO corrections to F2 calculated 
at twist-3 and taking M^v = 0, and we have compared it to the corresponding LO predic- 
tions. In this section we want to analyze how large are these corrections in comparison 
to other contributions: mass effects, higher-twists and corrections calculated at LO. 
These effects are calculated on the basis of App. [El and Ref. [2H]- 



6.2.1 F2 

The effect of including terms in LO twist-3 contribution to F2 corresponds to the 
change of the LO contribution obtained using the asymptotic DA for (— 3)-23 % when 
1 < < 10 GeV^. The change increases with and one obtains similar numbers for 
all three DAs. 

Let us now discuss higher-twist effects starting with twist-4. We note that for = 
the twist-4 contribution to F2 is (the whole contribution ( IE. lip . ( ]E.12p is proportional 



to Mjq). The ratio of LO twist-4 and twist-3 contributions to F2 is in the range (-5) 
to (-35) % for asymptotic DAs, 2 to (-0.7) % for BLW DAs, and 13 to 24% for CZ-like 
DAs, when 1 < < 10 GeV^. This behaviour is obviously very different for different 
DAs. In the case of the asymptotic DAs and the CZ-like DAs the absolute value of the 
ratio grows with Q^, while for BLW the ratio decreases, changes sign and then absolute 
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Table 7: The size of NLO twist-3 corrections to the proton form factor F2 (normahzed 
to LO twist-3 contribution calculated for Mj^ 7^ 0) in the range 1 < < 10 GeV^ and 
for three choices of scales: a) = = ^ = Q^, b) = Q^,fi% = /^r 1 = c) 



DAs 


I-iNLO,tw-3| 
^2 \Ml=0,a) 


7-ilNLO,tw-3| 
^2 \M%=Ofi) 


I-iiNLO,tw-3| 
-^2 Uf2^=0,c) 




T7LO,tw— 3 


T-iLOjtw— 3 
J^2. 


asy. 


62 - 57% 


71 - 36% 


58 - 51% 


BLW 


62 - 68% 


72 - 44% 


59 - 62% 


CZ-hke 


63 - 76% 


74 - 49% 


61 - 69% 



value increases. The role of the LO twist-4 contributions is thus small in the case of 
the results obtained using the BLW DAs, and more pronounced in the case of the other 
two investigated DAs (the BLW results seem to fall in some kind of minima). We have 
noticed the large sensitivity of these results on the choice of the parameters ff and A". 

Finally, we summarize our findings about the size of various contributions to F2 in 
Tables Eland [71 One can see that twist-3 contributions are dominant and positive. The x^- 
contributions are negative and the ratio of x^-contributions and LO twist-3 contributions 
does not change much for various DAs. The twist-5 contributions are more pronounced 
than twist-4 contributions and for both the ratio to twist-3 contribution is very sensitive 
to the shape of DAs. The twist-3 NLO corrections are positive and cca 60%. 

Hence, the twist-3 NLO corrections are indeed sizable and important. 

6.2.2 Fi 

In the next section we will proceed to the comparison of our results to the experimental 
data and for that we need the Fi contribution as well (we have at our disposal the 
experimental data for Gm and Ge)- Here we thus analyze the LO contributions to Fi. 

The effect of including terms in LO twist-3 contribution to Fi corresponds to 
the change of the LO contribution obtained using the asymptotic DA for 6-24 % when 
1 < < GeV^. The change decreases and then increases with (minimum at 
2 GeV^ and one obtains similar numbers for all three DAs. But in the case of Fi, 
twist-3 contribution is negative and small in comparison with twist-4 contribution. 

In contrast to -F2, twist-4 contribution to Fi is not proportional to Mjj (see App. 
IE|) . The effect of including terms in LO twist-4 to F2 corresponds to the change 
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Table 8: The size of various contributions to the proton form factor Fi (normahzed to 
LO twist-4 contribution calculated for Mjf ^ 0) in the range 1 < < 10 GeV^. 



DAs 


T-rLO.tw — 3 

fx 


T-rLO.tw — 5 

fx 


T-rLO.tw — 6 

f\ 


T-rLO. x'^-corr. 




T-iLO,tw-4 
^1 


T^LO,tw— 4 
^1 


j^LOjtw— 4 
^1 


asy. 


(-19) - (-7)% 


(-4) - (-5)% 


3 - 2% 


5 - 2% 


BLW 


(-25) - (-12)% 


(-2)% 


3-2% 


6-3% 


CZ-hke 


(-46) - (-105)% 


5 - 37% 


5 - 14% 


10 - 22% 



of the contribution obtained using the asymptotic DAs for 12-30 % when 1 < < 10 
GeV^. The numbers are similar for both asymptotic and BLW DAs, but smaller and 
negative for CZ-like DAs. The ratio of twist-3 and twist-4 LO contributions to F\ is 
(-19)-(-7)% for the asymptotic DAs, (-25)-(-12)% for the BLW DAs, and (-46)-(-105)% 
for the CZ-like DAs. Hence, apart from the results for CZ-like DAs at higher Q^, the 
twist-4 contributions is larger than the twist-3 contribution. 

Finally, to summarize our findings about the size of various contributions to F\ in 
Table [HI we state that the twist-4 contributions are dominant and positive. 

6.3 Comparison to experimental data 

Finally, let us compare our findings to experimental data. 

In Figs. [71 [Eland [9] we display the results for G m / {^^pG d) I-J'pGe/Gm and \/Q^ F2 / (KpFi) 
obtained using the asymptotic DAs (solid line), the BLW DAs (dash-dotted line), and 
the CZ-like DAs (dashed line). The DA parameters are given in Sec. 15.21 

For comparison, in Figs. [7^,[Hti' and[nK we present the LO predictions obtained on the 
basis of the results from Ref. [28], where the higher-twist contributions (up to twist-6) 
and x"^ correction were included and the value of Ai (16. 2p was used. The NLO predictions, 
i.e. the LO predictions obtained on the basis of the results from Ref. [28] plus NLO 
corrections for twist-3 (Mjy = case) calculated in this work (with sq and Mb as from 
Eq. fl6.ip . while fJ-ji = l^p = = ^s)' displayed in Figs. [7b, [Hb and [Hb- Here 
we investigate also the change of the results with the choice of Ai. The NLO results 
obtained using the default value of Ai (16. 2p are, as in Fig. [7tL,[Hti' and[9|L, denoted by thin 
lines, while thick lines denote the NLO results obtained employing the corrected value 

of Ai dnsi). 

In Fig. [7] we display the LCSR prediction for Gm normalized to fipGa. The displayed 
experimental data were obtained using Rosenbluth separation: A SLAC 1994 ^61j, ■ 
SLAC 1994 [62], ♦ JLab 2004 [63], * JLab 2005 [64j. The LO results [TJi favour the 
asymptotic and BLW DAs. The inclusion of NLO corrections (compare thin lines in 
Figs. [7K and [7b) raises the predictions. The change of Ai to the corrected value (16. 3p 
lowers the NLO results (thick lines) slightly. The NLO results seem to overshoot the 
data at lower Q^, while at higher again the asymptotic and BLW results seem to be 
closer to the data than the results obtained using the CZ-like DAs. 

In Fig. [HI we present the LCSR prediction for UpGe/Gm- We use the experimental 
data obtained using Rosenbluth separation ( ■ SLAC 1994 [62], ♦ JLab 2004 [63], 
JLab 2005 [M], A SLAC 1970 (small) [65J and Bonn 1971 (big) [HS] data revised in 
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(a) (b) 

Figure 7: LCSR prediction for the proton form factor Gm normalized to fipGu- The 
displayed results correspond to asymptotic DAs (solid line), BLW DAs (dash-dotted 
line), and CZ-like DAs (dashed line) - see Sec. 15.21 for the DA parameters, a) LO 
prediction on the basis of the results from Ref. [28] (i.e., the higher- twist contributions 
(up to twist-6) and correction included), b) LO prediction on the basis of the results 
from Ref. [2H] plus NLO correction for twist-3 {M^ = case) calculated in this work 
(choice of scales as in Fig. [6]). Thin lines: the default value of Ai (16.21) . Thick lines: the 
corrected value of Ai (16.31) . Experimental data obtained using Rosenbluth separation: A 
SLAG 1994 [61], ■ SLAG 1994 [62], # JLab 2004 [63], * JLab 2005 [M]- 




02468 10 02468 10 



[GeV^] q2 [GeV^] 

(a) (b) 

Figure 8: LGSR prediction for the proton form factors ratio ^pGe/Gm- The displayed 
results correspond to asymptotic DAs (solid line), and BLW DAs (dash-dotted line) - see 
Sec. 15.21 for the DA parameters, a) LO prediction on the basis of the results from Ref. 
[2H] (i-e., the higher-twist contributions up to twist-6 and correction included), b) 
LO prediction on the basis of the results from Ref. plus NLO correction (to F2) for 
twist-3 (Mtv = case) calculated in this work (choice of scales as in Fig. [6]). Thin lines: 
the default value of Ai (16.21) . Thick lines: the corrected value of Ai (16. 3p . Experimental 
data obtained using Rosenbluth separation: ■ SLAG 1994 [62], ♦ JLab 2004 [63], 
JLab 2005 [M], A SLAG 1970 (small) [65] and Bonn 1971 (big) [HB] data revised in [HZj. 
Experimental data obtained via polarization transfer: A JLab 2001 [68j, □ JLab 1999 
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(a) (b) 

Figure 9: LCSR prediction for the proton form factors ratio \/Q^F2/{KpFi). The dis- 
played results correspond to asymptotic DAs (solid line), and BLW DAs (dash-dotted 
line) - see Sec. 15.21 for the DA parameters, a) LO prediction on the basis of the results 
from Ref. (i.e., the higher- twist contributions up to twist-6 and x"^ correction in- 
cluded), b) LO prediction on the basis of the results from Ref. p8] plus NLO correction 
(to F2) for twist-3 (Mtv = case) calculated in this work (choice of scales as in Fig. 
Thin lines: the default value of Ai f l6.2p . Thick lines: the corrected value of Ai (16. 3p . 
Experimental data obtained using Rosenbluth separation: ■ SLAG 1994 [62j, ▲ SLAG 
1994 [61]. Experimental data obtained via polarization transfer: A and □ as in Ref. 
[28], Fig. 15 (M. Jones, private communication). 

|67] ) as well as more reliable experimental data obtained via polarization transfer (A 
JLab 2001 [UH], □ JLab 1999 The LO results displayed in Fig. [HJi show that 

while the results obtained using the GZ-like DAs are quite low and well beyond the 
data, and the results obtained the asymptotic DAs are on the high edge of the data, 
the BLW results seem to be in better agreement with the data, but one cannot really 
make some conclusive statements. The inclusion of NLO corrections (compare thin lines 
in Figs. [Hh' and [Hb) lowers the predictions significantly, while the change of Ai to the 
corrected value (16.31) raises the NLO results (thick lines) slightly. Note that the results 
obtained using the GZ-like DAs are ruled out and left out. One can see now the NLO 
results, especially the results obtained using the corrected value for Ai, are in quite 
good agreement with the data. The results obtained using the asymptotic DAs seem 
to describe well the experimental data obtained using the Rosenbluth separation, while 
the NLO results obtained using the BLW DAs seem to follow the slope of the preferred 
experimental data obtained via polarization transfer. 

In Fig. [9] we present the LGSR prediction for \/Q^ F2 / (KpFi) . We display the exper- 
imental data obtained using Rosenbluth separation (■ SLAG 1994 [62], ▲ SLAG 1994 
[51]) and preferred experimental data obtained via polarization transfer (A and □ as in 
Ref. [2H], Fig. 15 (M. Jones, private communication)). The LO results are displayed 
in Fig. [HK and while the results obtained using the asymptotic DAs are on the lower 
edge of the data, the BLW results seem to fall close to the data (at least for lower Q^). 
The inclusion of NLO corrections (compare thin lines in Figs. [9K and [9b) raises the 
predictions significantly, while the change of Ai to the corrected value (16. 3p lowers the 
NLO results (thick lines) slightly. As in the case of the ^Ge/Gm results displayed in 
Fig. [HI one can see that the NLO results, especially the results obtained using the cor- 
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rected value for Ai, are in good agreement with the data. Again, the results obtained 
using the asymptotic DAs seem to describe well the experimental data obtained using 
the Rosenbluth separation, while the NLO results obtained using the BLW DAs seem to 
follow the slope of the experimental data obtained via polarization transfer. 

In conclusion, the inclusion of NLO corrections calculated at twist-3 for = 
introduces significant changes in the LCSR predictions for Gm/{,1^pGd)i ^^pGe/Gm and 



\/Q^ F2 / (npFi) . It seems that NLO corrections, as well as the use of the corrected 

value for Ai (16.31) . bring the predictions for ii^Ge/Gm and \/Q^F2/{KpFi) in better 
agreement with the experimental data. For these quantities, the results obtained using 
the asymptotic DAs seem to describe well the experimental data obtained using the 
Rosenbluth separation, while the NLO results obtained using the BLW DAs seem to 
follow the slope of the experimental data obtained via polarization transfer. 

7 Summary and conclusions 

In this work the first attempt has been made to asses the size of NLO corrections to 
nucleon form factors. 

In LCSR approach dealing with nucleons is much more demanding than dealing with 
mesons, even at LO. For one, the number of contributing terms is rather large, the 
expressions are more involved, and the presence of three, instead of two, partons with 
corresponding momenta makes the calculation more complicated. All this is present 
at NLO also, with additional difficulties of one-loop calculation and larger number of 
contributing Feynman diagrams. 

In order to calculate the NLO corrections, we have started with the simple = 
(but Mn 7^ 0) approximation, i.e., the approximation corresponding to the first two terms 
in the expansion in nucleon mass (I2.44p . In that approximation only the leading twist, 
twist-3, and next-to-leading twist, twist-4, contributions appear. But to our surprise it 
turned out that the collinear divergences appearing in one-loop calculation do not cancel 
on the level of separate twist and that actually mixing appears which, without knowing 
the corresponding kernels, disables us in determing the finite contributions (see Sec. 14. 4p . 

Hence, we have strengthen our approximation and considered Mn = approximation 
(corresponds to the first term in (12.441) ) in which only twist-3 contributes and the evolu- 
tion kernels are known. We have shown the explicit cancellation of collinear, as well as, 
UV singularities (see Sec. 14.31) . The finite twist-3 NLO contributions to the correlation 
function are thus obtained in Mn = approximation and relevant invariant functions 
are listed in App. O 

We note that the observation of mixing of twist-3 and twist-4 NLO contributions is 
in nature similar to the observation given in Sec. 13.41 that the gauge invariant results 
are obtained not twist- by-twist but order-by-order in Mjy. The gauge condition is for 
Mjv = case satisfied both in LO and NLO order. For = and Mjy 7^ we have 
shown to LO that gauge condition is satisfied only when the sum of all contributing 
terms is taken into account, i.e., both twist-3 and twist-4 contributions. The additional 
condition is that the asymptotic forms of twist-3 DAs are used (no conditions, at least 
at this order, on twist-4 DAs). Hence, gauge invariance can be satisfied order by order 
in the expansion in M^r (I2.44p with possibly some additional conditions on the form of 




DAs. 
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To repeat, by switching-on the nucleoli mass, which is, of course necessary in order 
to determine higher-twists, we are at = and M^v 7^ stuck with mixing of the 
contributions corresponding to different twists. This "mixing" can be seen even at LO 
through the check of gauge invariance with respect to photon. If we take 7^ 0, there 
are no colhnear divergences whose cancellation we should take care of and no apparent 
"mixing". But the NLO expressions are more involved and there are additional open 
problems one should solve before attacking higher-order higher-twist calculations. For 
example, there is an open technical problem elaborated in Sec. 14.41 and connected to the 
calculation of the NLO contributions to second and third case defined in f l2.28p . The 
calculation of NLO corrections for 7^ case and thus NLO corrections to higher 
twists we postpone for some other time. 

In Sec. in] we have presented and analyzed our numerical results based on the cal- 
culation of NLO corrections in M^v = approximation, i.e., twist-3 NLO corrections in 
that approximation. Using loffe current in this approximation we are able to calculate 
only the corrections to F2 nucleon form factor. To make a full analysis and estimate the 
importance of NLO corrections, we have also included the LO results obtained beyond 
this approximation, i.e., leading twist and higher- twist results obtained in Ref. [2H]- We 
have considered here just the proton case. 

For Fi twist-4 LO contributions are dominant and positive, and there are no NLO 
corrections in = approximation. When one considers the size of various contribu- 
tions to F2, one realizes that twist-3 LO contributions are dominant and positive. The 
twist-5 LO contributions are more pronounced than twist-4 LO contributions and for 
both the ratio to twist-3 LO contribution is very sensitive to the shape of DAs. The 
^^-contributions are negative and the ratio of x^-contributions and LO twist-3 contribu- 
tions does not change much for various DAs. The twist-3 NLO corrections are positive 
and cca 60%. The NLO corrections to F2 calculated at twist-3 taking = are thus 
large. They vary for different DAs and depend on the choice of renormalization and 
factorization scales. In contrast to the dependence of F^^'-'/F^^, the dependence of the 
complete NLO prediction of F2 on the choice of renormalization and factorization scales 
is small. 

The inclusion of NLO corrections calculated at twist-3 for Mn = introduces signif- 
icant changes in the LCSR predictions for GM/ifJ'pGD), ^^pGe/Gm and \/Q^F2/{KpFi). 
It seems that NLO corrections, as well as the use of the corrected value for Ai (16. 3p . 
bring the predictions for ^pGe/Gm and \/Q^F2/{KpFi) in better agreement with the 
experimental data. For these quantities, the results obtained using the asymptotic DAs 
seem to describe well the experimental data obtained using the Rosenbluth separation, 
while the NLO results obtained using the BLW DAs seem to follow the slope of the 
preferred experimental data obtained via polarization transfer. 

Further analysis and inclusion of NLO corrections at higher-twists is needed to draw 
some more conclusive results. 
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A Feynman rules 

In calculating {—i T^) where is a correlator function (12.51) we use the standard Feynman 
rules for quark and gluon propagators and vertices, as well as for the quark-photon vertex. 
For loop integrals one has to introduce the usual integration over loop momenta^. 
The vertex corresponding to the interpolating nucleon current given by fl2.8p reads 

u-quark(a;, a)-u-quark(/3, 6)-d-quark(7, c)- loffe current : {C'^x)^^ (757'^) e"'"^ , 

(A.l) 

where all quark lines are going in the vertex and the order of u,u, and d-quarks with cor- 
responding Lorentz (a, /?, 7) and colour (a, 6, c) indices is counterclockwise. Furthermore, 

(C^7a)^ = C^7A, i.e., (C7a)„, = {C^>)p^. 

The (incoming) nucleon "projector" corresponding to (12. lip and the first case of 
^TM is given by 

u-quark(MiP, a, a)-u-quark('U2-P, &)-d-quark('U3P, 7, c): 

- j VuF^HuuU^u,) X% F« — , (A.2) 

where all quark lines are outgoing from the nucleon blob and the order is clockwise. The 
trivial identity xf^ = (X^^)) together with (^J^ is useful m some cases. 

The typical contribution obtained using the general Lorentz decomposition (12. lip 
and loffe current (12. 8p . i.e., Feynman rules ( 1A.2P and ( lA.ip . respectively, has two parts. 



For the d — d quark line, by going, following the standard rule, in the opposite direction 
of the fermion line, one obtains the product of 7 matrices with the nucleon spinor. The 
u — u lines close the trace and obviously, in writing it down, one goes opposite to the 
direction of the one quark line and along the other one. The latter case corresponds to 
{l^lll^l2 ■ ■ ■ Ifin)'^ (where 7^1 . . . 7^^ is the order of 7 matrices opposite to the direction of 
the fermion line) and one then makes use of 

{1,11,2 ...i,y = {-irc^,n . . . i,,c-' , (A.3) 

i.e. when going in the direction of fermion line one puts the 7 matrices on the that line 
between C and C~^. 

Finally, let us mention that the usual relations for S\]{Nc) algebra should be employed 
in calculating the colour factors. Since we are dealing here with nucleon described by 
three quarks we are actually already assuming Nc = 3 and only this choice leads to 
gauge invariant results. Obviously, 



^^We use dimensional regularization in 13 = 4 — 2e dimensions and for integral measure we choose 
IJ?'^ J d^l(2-K)^ - see, for example, App. C in Ref. [5^ for some comments on this choice, i.e., on 
introduction of scale in Feynman integrals. 
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B 75 ambiguity in dimensional regularization 



When using dimensional regularization, one runs into trouble with quantities that have 
the well-defined properties only in D = 4 space-time dimensions, that is, with the Levi- 
Civita tensor e^^^xK, which is a genuine 4 dimensional object, and consequently with 
the pseudoscalar 75 Dirac matrix. Let us mention that the appearance and mixing 
with evanescent operators [TO], as well as, the definition of Fierz transformation in D 
dimensions are also connected to this problem. We shall handle it similarly to Ref. [52] 
with some additional finesse concerning Chisholm identity (see Sec. 14.21) . Below we 
explain the general features of the ambiguities that we encounter in our calculation. In 
order to resolve these one should generally use some other input like the knowledge of the 
quantity that does not suffer from ambiguities, condition of cancellation of singularities, 
condition of preservation of gauge invariance. Ward identities etc. 

B.l General remarks — trace ambiguity 

The generalization of the 75 matrix in D dimensions represents a problem, since it is not 
possible to simultaneously retain its anticommuting and trace properties. In practice, 
the ambiguity arises when evaluating a trace containing a 75 and pairs of contracted 7 
matrices and/or pairs of Dirac slashed loop momenta . To deal with a 75 matrix, several 
possible schemes have been proposed in the literature. 

In the so-called naive-75 scheme ^3], the anticommutation property of 75 is retained, 
while the cyclicity of the trace is abandoned. The traces obtained by cyclic permutation 
of the matrices differ hj D — 4. Consequently, if the trace is multiplied by a pole in D — A, 
there appears a finite ambiguity in the result. An alternative scheme has been proposed 
in the original paper on the dimensional regularization by 't Hooft and Veltman [55] , 
and further systematized by Breitenlohner and Maison [32] • In this scheme, to which we 
refer as HV scheme, the anticommutativity of 75 is abandoned. In contrast to naive-75 
scheme, this scheme is claimed to be mathematically consistent but still not without 
drawbacks. Namely, this prescription for 75 violates the Ward identities and introduces 
"spurious" anomalies which violate chiral symmetry. To restore the Ward identities, 
finite counterterms should be added order by order in perturbation theory [7T]. In this 
scheme, the cyclicity of the trace is retained. 

If a trace contains an even number of 75 matrices, then the property 75 = 1 can be 
used to eliminate 75's from the trace, and the Ward identities are preserved if the naive- 
75 scheme is used [S^ (the cyclicity of the trace is restored and the corresponding results 
are unambiguous). On the other hand, in the HV scheme the "spurious" anomalies can 
occur owing to the non- anticommuting property of 75. As for the traces containing an 
odd number of 75 matrices, we are left with the above mentioned ambiguities in the 
results. 

For details, we refer to App. A in Ref. [52] . 

B.2 General remarks — Chisholm identity 

Additionally, the Chisholm identity that we need in our calculation 
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is strictly speaking valid only in D = 4 dimensions. The modification for HV scheme 
can be found in the literature (see, for example, Tracer [72] manual) but we have not 
been able to find any recipe for the naive-75 scheme. 

Our analysis of the problem has shown that when applying the Chisholm identity on 
expressions of the form 

Yl.Ye>..ai3 (B.2) 

different results appear in dependence of whether one first contracts the Levi-Civita 
tensor with the 7 matrix on the left or the right side of the non-contracted matrix (7^; 
in this case). The difference is again, as expected, proportional to D — A. For example, 

7'^7"7k7'^7^^m-/3 (B-3a) 

has two sets of results 

(7^^^) 7"7.7'7^ = 1^ ire^uap) 7.7^7^ = -^{D - 4){D - 2){D - 1)7,75 , (B.3b) 
while 

7'^7"7. (Ye^uaf^) 7^ = 7'^7"7«7' {l^e,,^^) = i{D - 4)(D - 2){D - 1)7,75 . (B.3c) 

So, when using the Chisholm identity in its form as in Z) = 4 dimensions (which should 
be in agreement with the "philosophy" of naive-75 scheme) we again, as in the case of 
trace ambiguity, encounter the ambiguity proportional to D — 4, which, when multiplied 
by pole in D — 4, possibly leads to finite ambiguity of the results. 

C NLO results for Mat = case 
C.l Formalism and notation 

In M^r = case we cannot asses the A contribution from Eq. (12.341) . while the corre- 
sponding B contribution we list in this section. 

The B function can be written as a convolution in terms of Vi and Ai nucleon DAs 

^Mjv=o(Q^) P'^) 

= ^^3,yl,MJv=o({a;fc},Q^P'^;/i|) ® Vl({xfc};/i|) 

+TB,Ai,A%=o(K},g',P";/i|) ® Ai({xfc},/i|), (C.l) 

or in terms of twist-3 nucleon DA $3 = Vi — Ai 
Bmi^=o{Q'^, P'^) 

= Ts,^,MN=oi{xk},Q^P'^;f^l) ® $3(K};/x|). (C.2) 

Here q = —Q^ is a photon virtuality and P' = P — q, while P is incoming nucleon 
momentum. The factorization scale is denoted by fi'jp, and by {xk} the quark momentum 
fractions xi, X2 and 0:3 are denoted. Note that xi and X2 denote the momentum fractions 
of u-quarks, while X3 correspond to the momentum fraction of the d-quark. A stated in 
Eq. (12.251) . Vi is symmetric and Ai antisymmetric under xi X2 exchange. 
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It is convenient to introduce a dimensionless quantity 

and we will from now on expres^ the B function using W instead of P'^. 
Generally, we write 

with F^*^ denoting nucleon DAs and ^'^^^(i) j^^^^q the corresponding "hard-scattering" 
part. The nucleon DAs are intrinsically non-perturbative quantities but their evolution 
to scale can be calculated perturbatively. Nevertheless, we take into account only 
the LO evolution or neglect the evolution completely. The "hard-scattering" part is 
calculated perturbatively and in this work the calculation to NLO is performed. Hence, 
following SecO (see Eqs. and IK^ . and Eqs. (KM and (OTj) ) we can write 

the expansion of T^^F(i)^MN=o 

{{xk},W-fil/Q') 

(C.5) 

where 

and /ifj and //^ ^ scales denote the coupling constant and loffe current renormalization 
scales (which are in practice often taken the same, and even same to the factorization 
scale fip, but are in principle independent). 

Note that all order result for Tg p{i) Mm=o would not depend on the choice of the 
renormalization scale, but the truncation of the series to any finite order (in or case 
NLO) introduces the residual dependence. This dependence would be stabilized by 
inclusion of higher-orders (a", n > 2). One is left also with the residual dependence of 
B on the factorization scale (see Ref. [73] for details on that point). 

In the following we summarize the LO and NLO results for Tq pn) Mn=o- These are 
proportional either to e„ or being u and (i-quark charges (depending on where the 
photon coupled), respectively. Remember that we are displaying here the proton case. 



Although, the functional dependence on W is different from the one on P'^ we retain the same 
nomenclature, i.e., from now on we use B'{Q^, W) = B{Q^, W). 
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while for the neutron e„ and have to be exchanged. Hence, we burden our notation 
with one more index 



T, 



(...]= T" '^"^ (■ 



+ T'" 

I YD 



(■■■)• 



(C.7) 



In order to simphfy the expressions and write them in a form most suitable for further 
calculation, we introduce the following functions 



97 {Xi, Xj 
dsiXi, Xj 
ggi^Xi, Xj 

gioixi,Xj 
gn{xi,Xj 

9l2{Xi,Xj 



W) 
W) 
W) 
W) 
W) 
W) 
W) 
W) 
W) 
W) 
W) 
W) 
W) 



1 



{xiW -l+irj)' 

ln(l — XiW — if]) 
(l-XiW-iri) ' 

ln^(l - XiW - if]) 



{1-XiW - 


-i7]) 


ln(l - XiW - 


iri) 


{W + IT]) 




ln(l - XiW - 


ir]) 


{W + ir])'^ 




\n^{l - XiW - 




{W + 17]) 




\n^{l-x,W - 


-ir]) 



{W + tr])^ 

ln(l — {xi + Xj)W — ir]) 
(1 - XiW - ir]) ' 

ln^(l - {xi + Xj)W - ir]) 
(1 - XiW - ir]) 

ln(l — (xj + Xj)W — ir]) 





{W + ir]) 




ln(l - 


- {xi + Xj)W - 


-ir]) 




iW + lT])'^ 






- {xi + Xj)W 


- ^V) 




{W + ir]) 






- {xi + Xj)W 


-ir]) 



{W + ir])^ 



(C.8) 



Note that we have kept ir] terms (?7 > and r] <^) coming from the Feynman diagram 
calculation (quark and gluon propagators), which will enable the correct determination 
of imaginary parts necessary for LCSR in Sec. O 



52 



C.2 Complete list of results 

In previously introduced notation the LO contributions to Tq y^^Mn=o Tis^Ai,Mn=o 
read: 

7S,%V=o({^J,W^) = 0, (C.9) 

and 

TBy,M.=o(M,W) = e^[go{xuW)+goix2,W)], 

TBAuMr.=oiM.W) = eu[-goix,,W)+goix,,W)]. (C.IO) 

Furthermore, the NLO contributions proportional to ln(yu|j ^/Q^) take also the simple 
for 



and 



TZTZt'oiM^W) = 0, (C.ll) 

r?yX=o(K}.^) = 2e^ [go{x,,W) + goix2,W)] , 
^SXlolW^W^) = 2eu[-goixuW)+goix,,W)]. (C.12) 

Next we turn to NLO contributions proportional to ln(/i|,/(5^) which originate from the 
factorization of collinear divergence^ and read 

Tsy?:MfMM,W) = p,{Qgoix,,W) + 8g,ixs.W) 

4X1X2 - X1X3 - X2X3 X1+X2 , 



-2 ^7(^3, xi, ly) - 2 ^7(^3, X2, W) 

^ 9%{^Z, Xx,W)- — gg{x3, X2, W) 

Xi X2 

gio{x3,Xi,W) 7 gioix3,X2,W) 



a;i(a;i+a;3) X2\X2^X'i) 
TtS:U{^.}.W) = UJ^^9s{xs,W)-^-^g,{x„W) 

O ( XlX0X3 



X1X2 X1X2X3 

+ — gg{x3, Xi, W) - — 5(9(^3, X2, W) 
Xi X2 



^^They are as expected proportional to LO and loffe current renormalization factor 2 - as renormal- 
ization of UV divergences demanded (see Sec[ 



15 



For the details we again refer to Sec 
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9W{X3, Xi, W) H r gw{x3, X2, W) 



and 



Xi{xi+X3) ' ' X2{X2 + X3) 



+4 g^{x^,W)+Ag,{x2,W) 

Xi - 2X3 I Tj.x , X2 - 2X3 I Tj.x 
X\X3 X2X3 

xia;2 XxX2 

- gi{xx.,x2.,w)- 5'7(x2,xi,iy) 

- 5(7(2:1, a;3,iy) - 5(7(2:2, xs, PF) 

— ^ S-gla^i, 2^3, Vl^) - — 5^9(3:^2, X3, W) 
X3 X3 

g\^{Xx,X2,W) 



X1X2 



-4 gi(x,,W)+Ag,(x2,W) 



Xi — 2X3 X2 — 2X3 

gsixi, W) + 53(2:2, W) 



X1X3 X2X3 
^ g4{xi,W) + ^ g4{x2,W) 



(C.13) 



a;ia:;2 X1X2 
+ g7{xuX2,W) - g7{x2,Xi,W) 

+ 57(2:1,2:3,1^)- 57(2:2, 2:3, PF) 

+— gs{xi, X3, W)- — gg{x2, X3, W) 
X3 X3 

+ ""r":' , gio{xi,x2,w)] . (C.14) 

2:12:2 (2:1 + J 

Finally we give the lengthy expressions for "finite" NLO contributions 
TSy?:M7=oiM:W) = ^e,{l8go{x3,W) + 12g,{x3,W)-8g2{x3,W) 

12 , T^^s 5(,Xi+2:2) , 

— g3(x3, w) + ^ ^ ^4 2:3, w) 

X3 XIX2X3 

4XiX2 - X2X3 - X1X3 , , 2:1 + X2 , 

+ 55(2:3, W) gG{x3, W) 

XIX2X3 XIX2X3 
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-3 g-jix^, xi, ly) - 3 g-j{x^, X2, W) 

+2 58(a;3, xi, W) + 2 gs{x3, X2, W) 

5 5 

gio(xi,X3, W) ■ gio(x2, X3, W) 



Xi{xi + X3) ' ' a;2(a;2 + a;3) 

+— guixi, X3, W) + — gu{x2, X3, W) 

Xi X2 

guixi, X3, W) H gi2{x2, X3, W) \ , 



Xi{Xi+X3) X2[X2+X3) 



(C.15a) 



X2 — Xi Xi — X2 I . 

X\X2 X\X2X3 

-3 g'j{x3, xi, ly) + 3 g'j{x3, X2, W) 

5 5 

5fio(xi, X3, W) H ■ giQ{x2, X3, W) 



Xi{xi+X3) X2{X2+X3) 

— ^ gi\{xi,X3,W) + — gil{x2,X3,W) 

Xi X2 

1 1 

512(2^1, X3,W) ■ ^ 512(2^2, 2^3, 



Xi{xi+X3) ' ' X2{X2+X3) 

(C.15b) 



and 

nNLO,fin,e 



2 C 19 19 
= ^(^u{-^ go{xi,W)- — gQ{x2,W) 

-6 gi{xi, W)-Q gi{x2, W) + 4: g2{xi, W) + 4: ^2(2:2, W) 

2xiX2 + 6X2X3 - 3xiX3 2xiX2 + 6X1X3 - 3x2X3 

H 53(2^1, W) + g3{x2, W) 

X1X2X3 X1X2X3 

4 4 

giixi,W) g4ix2,W) 



X1X2 X1X2 

, Xi - 2X3 X2 - 2X3 

H gbixi, W) H g5{x2, W) 

X1X3 X2X3 

+^ ge{xi, W) + ge{x2, W) 
X1X2 X1X2 

+3 57(2^1, X2, W) + 3 g7{x2, xi, W) 

- g8{xi,X2,W) - g8{x2,xi,W) - g8{xi,X3,W) - g8{x2,X3,W) 
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_^ 3{xi + X2) ^^^^^^ - A g^(^xi, Xs, W) - — g9{x2, X3, W) 

X1X2 X3 X3 

4 

H gw{xi,X2,W) 

X1X2 

— - gii{xi,X3,W) - — gii{x2,X3,W) 
X3 X3 

^ gi2{xi,X2,W) \ , (C.16a) 



X1X2 



2 f 23 23 

= :^eu<Y9o{xi:W)- — go{x2,W) 

+6 gi{xi, W)-6 gi{x2, W) - A g2{xi, W) + A g2{x2, W) 

, 4xiX2 - 6X2X3 - 70:1X3 , ^j^, , -4xiX2 + 6X1X3 + 7x2X3 , ^j^, 

+ ^3(2:1, W) + g3{x2, W) 

X1X2X3 X1X2X3 

2(x2 - 5x3) 2(xi - 5x3) 

^4(3:^1, W) + g4{x2, W) 



X1X2X3 ^1X2X3 

2X3 - Xi X2 - 2X3 

H ^5 a^i, W) H W^) 

X1X3 X2X3 

— ge{xi, W) + ge{x2, W) 

X1X2 X1X2 

-3 g7{xi,X2, W) + 3 g-j{x2, xi, W) 

+ gs{xi,X2,W) - gs{x2,Xx,W)+ gs{xi,X3,W) - g8{x2,X3,W) 

j^JiM ^ ^^^^^^ _ A ^g^a;^^ 2:3^ + A ^g(^a;2, X3, 

X1X2 X3 X3 

10(xi - X2) , ^j... 

gio{xi,X2,W) 



XiX2(xi + X2) 

2 2 

H 7 ^ r giQ{.Xi,X3, W) ■ 5^10 (a;2, 2:3, W) 

X3(Xi+X3) X3(X2 + X3) 

+ — gi\{xi,X3, W)- — 5fii(x2, X3, W) 

X3 X3 

^^i-X2 g^^^x^^x^^w)] . (C.16b) 



XiX2(xi + X2) 

The contributions to $ M;v=o which, as can be seen, due to symmetry properties of 
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contributions Tg,vi,Miv=o and Tq^Ai,Mn=ol3, correspond exactly to 



B,<S>,Mn=0 — Tt3yi,MN=0 — -iB,Ai,M]v=0 



(C.17) 



C.3 Summary for proton case 

Finally in this section we list the expressions that we actually use in our numerical 
calculations of the proton form factors (see Eqs. ( 1C.2IIC.8I) for corresponding definitions). 
We have already made use of e„ = 2/3 and = —1/3. 
For Mat = twist-3 LO contribution reads: 



NLO contribution proportional to ln(yu|, ^^/Q^) is given by: 



(C.18) 



'^go{xi,W) + ^go{x3,W) 



NLO contribution proportional to ln{fijp / Q"^) reads: 



T-iNLO,IR /r 1 Ti/\ 

= 4 



4 2 16 8 

- goixi, W)-- goix3, - — gi{xi, W) - - gi{x^, W) 



4(2a;3 - xi) 2(2xi - X3) 



9x1X3 



9x1X3 



9XiX2 



4 2 

g^ixi, W) - ^4(3:3, W) 



9x2X3 



(C.19) 



4 4 2 2 

+ g g7{xi,X2, ^) + g 97{xi, X3, W) + - gj^xs, xi, W) + - gr^xa, X2, W) 

, 2(2x1 + X3) . .... 
+ — ^9(3^1, a;3, W) 



9x1X3 
4 



9X2(Xi + X2^ 



■gw{Xl,X2,W) + 



9X2 (X2 + 2:3) 



gio{x2,xs,W) 



(C.20) 



"Finite" NLO contribution reads: 

{{xk},W) 

= 4 



7 2 

- 3 go{xi, + - go{x2, W) - go{x3, W) 



giixi, W)-^ giixs, W) + ^ g2{xi, W) + ^ (72(^3, W) 

+ + 93{xu W) + l(^-^) g,{x2, W^) + gsixs, W) 

9 \Xi X2 X3 / 9 \X3 Xi / 3x3 



"'^^The former is symmetric under xi <-> X2 exchange and the latter antisymmetric, analogous to 
symmetry properties of Vi and Ai . 
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2{x2 - 7x3) 2(xi - 3x3) 5 

+ ^^^ ^^4(3:1, W) - ^^4(2:2, W) - ^4(3^3, W) 

9x1X2X3 9x1X2X3 9x2X3 

2(xi — 2x3) , (x3 — 2xi) , 

+ ^7^ -95{xi, W) + ^ '-g,{x3, W) 

9x1X3 9x1X3 



2 2 1 

-- gfi{xi,X2, W^) - g gfi{xi,X3, 

2 



5 



9 \xi X2, 



-^g7{x3,X2,W) 

2 2 

I 5(9(xi, X2, W") + — 5f9(xi, X3, W) - —gg{x2, X3, W) 

ilAxi — 6x2) T,,. 2 

+ 7^ 7 \ ^gw{Xi, X2, W) - — — ^ -^io(Xi, X3, W) 

9XiX2(Xi + X2) 9X3(Xi + X3) 

, (2x2 + 5x3) (2xi + X3) 
+ 7^ 7 \ ^9w[x2, X3, VT) X3, W) 

9X2X3(X2+X3) 9X1X3 

2 1 

7 r5'i2(a;i, X2, W^) 7 r5'i2(a;25 

9X2(Xi+X2) 9X2(X2+X3) 



.X3,W) 



(C.21) 



D Imaginary parts of selected functions 

In this section we list the imaginary parts of selected functions that appear in our cal- 
culation. 

We start with the well-known result for the logarithmic function 

ln(x — Xo ± irf) = ln(|x — Xo|) ± in 0(xo — x) . (D.l) 

Here x, xq, and t] are real, and 77 > 0, -C. All other results that we list can be derived 
from (iDlll . 

It follows trivially that 

ln^(x — Xo ± = [ln^(|x — Xo|) — 7r^6(xo — x)] ± 2i7i ln(|x — Xo|)6(xo — x) , (D.2) 

and higher exponents ln"( X — Xq ± irj), n > 2 can be obtained similarly. 
The other well-known result 

-— = V =F ^vrdfx — Xo) 

X — Xq ± «?7 X — Xo 

= [9(x — Xo) + 0(xo — x)] =F — a^o) (D-3) 

X — Xo 

can be obtained by taking a derivative of Eq. fID.ip with respect to x: 

= — ln(x — Xo ± irj) . 

X — Xq ± IT] ax 
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Furthermore, since 



1 



(x — xo ± ir])" (n— 1)! dx" ^x — xo±irj 
it is easy to see that for n > 2 
1 



(x — xo ± ^77)" (x — Xo) 
Finally, by making use of 



^ [9(x - Xo) + e(xo - x)] T iTT 1„ i <5("'^)(x - Xo) ■ (D.4) 



-— ln^(X±^r^) 



(n- 1)! 



ln^(X±ir7) 1 d 



X±iri 2dX 
where X = x — xq or X = xq — x, we get 

9(xo - x) 



X±iri 3dX 



-— ln^(X±^r^) 



and 



1 ^ ln(x — Xq ± tri) 
— Im 

vr X ~ Xq ± ir] 



1 ^ ln(xo — X ± zr?) 
— Im 

TT Xo — X ± ir] 



± 



± 



± 



± 



X — Xq 



— 6{x — Xo) ln(xo — x) 



J 6(xo — x) 
\ X - Xo J + 



5(x — Xo) ln(|a — Xo|) 



Q{x - Xo) 
Xo — X 

6(x — Xo) 
Xo - X J ^ 



— 6{x — Xo) ln(x — Xo) 

- 6{x - Xo) ln(|xo - b\) 



1 In (x — Xq ± iri) 
— Im 

71 X — Xq ± ir] 



± 
± 



6(xo — x) 



2 ln(xo — x) 

X — Xo 



+ 6{x — Xo) 



IT 



ln^(xo — x) 



e(xo-x)^^^^^^^^|> +5(x-Xo) 



X — Xo J ^ 



- ln^(|xo - a\) 



(D.5a) 



(D.5b) 



(D.6a) 



1 ^ ln^(2;o — X ±iT]) 
— Im 



TT 



Xq — X ± if] 



± 
± 



0(x — Xo) 



2 ln(x — Xo) 



e(x - Xo) 



Xo — X 

2 ln(x — Xo) 



+ 5{x - Xo) 



— ln^(x — Xo) 



+ 6{x - Xo) 



TT 



y-ln2(|xo-6|) 



(D.6b) 



Note that right hand-side of the first lines in ( ID.SP and ( ]D.6I) consists of two terms which 
separately "blow up" for x ^ xo- The sums are finite and in the following lines, using 
the usual {}+ prescription 



{F(x,xo)}+ = -F(x,xo) - (5(x - Xo) / dzF{z,xo) 



(D.7) 
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we express them as manifestly finite sum of two terms. It is easy to see that 




X — Xq J _|_ X — Xq 

6(x — Xq) 



— 6{xo — x) ln(xo — Xq) + 6{x — Xq) ln(|a — xo|) 

— S{x — Xq) ln(xo — Xq) + S{x — Xq) ln(|a;o — b\) , 



Xq — X 

(D.8) 

and similarly for < B(x — xo)^^^^^^r^ r \ ©(^^o ^ x) \ . Note that in this 

section we take integration variable of the next step of the calculation and the 

variable with respect to which the {}+ prescription has been defined^. 



E LO results and LCSRs for 7^ case 

In this section we present LO twist-3 and twist-4 contributions to A and B functions 
fl2.34l) based on results listed in Tables [21 E] and H] {Cp^jJ^ and Cp^^^ coefficients, 
respectively) which were calculated for Mn ^ 0. Furthermore, using (12.431) we determine 
the LCSR contributions to form factors Fi(Q^) and F2((5^). 
Analogously to Eq. fIC.SI) . we define 

^/o(x„iy,M^/Q2) ^ _1 _ (E^^) 

[XiW - 1 - Xi(l - Xi) + IT]) 

Note that g^ixi^W) introduced in fIC.Sp corresponds to g^ixi^W,^), i.e, in other words, 
in (lE.ll) we introduce the generalization to 7^ 0. In this calculation we will also 
encounter 

gl{x,,W,MllQ'') = . (E.2) 

{XiW - 1 - ^Xi{l - Xi) + 17])'' 

The twist-3 contributions and the twist-4 contributions corresponding to nucleon DAs 
V3 and A3 can be expressed in a form of a convolution: 

^LO(Q^ly,M^;/x^) = ^ E T^%i{xk},W,M'^/Q')0F^^i{xk};f^l), 

SLO(g^ W, M^; /i|) = ^ T^L° „ ({x,}, W, Ml/Q') F« ({x,}; /i|) , (E.3) 

^ i 

where F« e {V^i, Ai, I/3, A3}. 

The coefficients Cp°^ from Table [5] determine the "hard-scattering" twist-3 con- 
tributions to A function 



TA%iM,W,Ml,/Q^) = -2eu[xigo{xi,W,Ml,/Q^)+X2go{x2,W,M^^/Q' 
T^%,i{xk},W,Ml,/Q') = 0, (E.4) 



"'^'''Some similar results but not for general {}+ prescription (a = 0) can be found in App. A of Ref. 
S]- Note the typo in Ref. [32]: 6-function is not written inside the {}+-prescription. 
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while the coefficients C'p^^^' ' determine the "hard-scattering" twist-3 contributions to 
B function: 



T^%{{xk}.W,MllQ^) = eu[-g^{x^,W,MllQ^)+g,{x2,W,MllQ^)\ . 

Similarly from Table [3], the "hard-scattering" twist-4 contributions to A function are 
given by 

T\%,{{xk},W,Ml/Q'') = ?,eu[xig,{xi,W,Ml/Q'')+X2go{x2,W,Ml/Q'')\ 

+2edX^g^{x^,W,Ml/Q^), 
T\'^A,{{xk}.W,Ml/Q^) = e4-xi (7o(xi, PF, M^/g2) + (7o(x2, VT, M^/Q2)] (E.6) 
while the contributions to B function read 

Th%{{xk}.W,Ml/Q^) = 0, 

T^%{{x,},W,Ml/Q^) = 0. (E.7) 
The contributions listed in Table H] contribute to functions A and B according to 

where F123 ^ {^123)^123} ^ corresponding definitions fl3.15H3.l5I) . 
The "hard-scattering" contributions to A then read 

T\%,,{X2, W, Ml/Q^) = eu [gl{x2. W, Ml/Q^) - go{x2, W, M^/Q^)] , 

T^y,Jxs, W, Ml/Q^) = 2e, [g',{xs, W, Ml/Q^) + g,{x,, W, M^/Q^)] , (E.9) 

and 

T'^A.^^^Xr.W^MllQ^) = -e4glix^,W,Ml,/Q')+goixuW,M'M/Q')] , 

T!fA,J^2, W, Ml/Q'') = e„ [gl{x2, W, Ml/Q'') + g,{x2, W, M^/Q^)] , 

T\X->^^^^Mn/Q^) = 0. (E.IO) 
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The "hard-scattering" contributions to B are given by 



/Vf2 



T^%Jx,,W,M'^/Q') = 2e,-^xsgl{x,,W,M'^/Q'), (E.ll) 



and 



/bf2 

TkXrM^W,M'^/Q') = -e^-^x,gl{x,,W,Ml,/Q'), 

TkXjx2,W,Ml,/Q') = e^-^x,gl{x2,W,M'^/Q'), 

T^Xjx,,W,M'^/Q') = 0. (E.12) 

For Mjv = there are no contributions to A since in the decomposition of the 
correlation function f l2.34p A is multiphed by Mn- Moreover, in the hmit only 

the twist-3 contributions to B, given in ( 1E.5I) . "survive" and take the form ( IC.QIIC.lOj) . 
By taking M/v 7^ but = (corresponds to the first two terms in the expansion 
(12.441) ). one is left with the same twist-3 contribution to B ( (lE.Sp with = 0, i.e., 



f lClMcnol) ) and twist-3 and twist-4 contributions to A: (El), (EE]), (lK9]) and flKTOl) 
taken with = 0. Note that in comparison to twist-3, the twist-4 contribution to 
B are suppressed by Mj^/Q^ factor (see flKTl) . flKTTl) . flKT2l) ). This is not the case for 
A (see (IE.4p . flE.6p . flE.9p . flE.lOp ) where there is no such additional suppression factor 



between twist-3 and twist-4 contributions. For Mjf 7^ the higher twists also contribute 
(twist-5, . . . ), which we do not consider here but just refer to the results presented in, 
for example, Ref. [28] (App. A). From these results one can see that both twist-4 and 
twist-5 contributions to B are suppressed by M^/Q"^ in comparison to twist-3. For A, 
twist-5 and twist-6 contributions are suppressed by / in comparison to twist-3 and 
twist-4 contributions (remember that there is an additional factor M^r infront of A). 

According to Eq. (12.431) and analogously to (15.51) . we can now formulate the "rules" 
for separate terms contributing to T^^ and T^^ leading to separate terms contributing 
to Fi and F2, respectively. The contributions (1E.4I - IETTI) can be conveniently expressed 
as a sum of the terms of general form 

go{{xk},W,M'^/Q^) /(K}), (E.13) 
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and they then contribute as 



^ ' Vx dw e-(-i)QVA^|+M^/Af: 



AiTT 



X lm[go{{xu},w,Ml/Q^)\ /({x,}) 



XQ Jo 



X 



X 



Here F» e {Vi, Ai, V3, A3} and 

V(Q2 + So-M2,)2+4M2,Q2 _ (Q2 + _ M 



Xo 



(E.14) 



(E.15) 



Note that hm^j2^^o^o = Q'^/iQ'^ + -So), i-e., one recovers the lower hmit from Eq. (15.51) . 
The contributions ( lE.9IIE.12l) consist of the terms of the form 



go{x,,W,Ml,/Q^) fix,) and g^,{x,,W, M^^/Q^) /(x, 
Analogously to ( IE. 141) . the former contribute as 



(E.16) 



LO,{Vi23,Ai23}/^2. ,,2 



,LO,W23,Ai23}^g2.^2^) 



1 r('o+Q^)/Q'' 



dxk 



X 



Im [5(o(xfc,w,M^/g^)] f{xk) Fi!^3{xk;id%) 



- t dxk e-(i-^fe)QV(^'fcAf|)+xfeM2,/M; 



Ai 



1 



(E.17) 
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The terms with take shghtly more comphcated form 



2F[ 



LO,{Vl23,Al23}/^2. ,,2 



1 



"1 r{so+Q^)/Q'^ 
dxk 

Aivr Jo Ji 



X 

1 



Im [gl{xk,w,M^/Q^)] f{xk) FS^{xk; nl) 



Q 



-(so-Ml)/Ml 



g2 + xgM 



/(a^o) ^123 (^^fe = a;o;/w|) 



AT 



Ml J., 



Xi 



(E.18) 



where F«3G{ni, 433}- 

We note that these resuhs are in agreement with the somewhat differently derived 

expressions from Ref. [28 j. Eqs. (A.15-A.18. [28J). We refer to that paper for higher-twist 

contributions. 
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